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ABSTRACT 


The  application  of  Direct  Bit  Detection  methods  are  analyzed  and  evaluated  in 
digital  communication  systems  employing  Differential  Phase  Shift  Keyed  (DPSK) 
modulation.  Assuming  an  additive  white  Gaussian  noise  interference  model,  and  direct 
bit  detection  methods  used  in  the  receiver,  8-DPSK  and  16-DPSK  communication 
systems  are  considered  and  their  performances  evaluated  in  terms  of  the  delivered  bit 
error  rate  as  a  function  of  signal-to-noise  ratio.  The  advantages  and  disadvantages  of 
Direct  Bit  Detection  Receiver  (DBDR)  systems  used  in  conjuction  with  differential 
phase  encoding  are  d  ermined,  with  specific  application  to  8-DPSK  and  16-DPSK 
modulated  signals.  The  effect  of  improper  receiver  carrier  phase  synchronization  is 
considered  in  detail,  and  resulting  performance  degradations  are  evaluated.  Numerical 
results  show  that  for  receiver  phase  errors  of  more  than  a  few  degrees,  severe 
performance  degradations  result  for  both  DPSK  modulation  schemes,  unless  a 
complete  phase  reversal  (i.e.,180°)  takes  place. 
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I.  INTRODUCTION 


M-ary  differential  phase  shift  keyed  (M-DPSK)  modulation  is  a  widely  used 
technique  in  digital  communication  applications.  Bandwidth  efficiency,  relatively  good 
noise  immunity,  constant  signal  envelope,  and  simplicity  of  implementations,  make  this 
scheme  particularly  attractive  for  use  over  satellite,  terrestrial  radio  and  voiceband 
telephone  channels.  While  system  analyses  pertaining  to  the  performance  of  M-DPSK 
receivers  abound  in  the  literature,  treatment  is  usually  restricted  to  the  case  where 
signal  transmission  takes  place  over  an  additive  white  Gaussian  noise  channel.  [Ref.  lj 

Since  the  digital  information  transmission  is  based  on  carriers  modulated  by  the 
symbol  waveforms,  the  symbol  error  probabilities  which  quantify  receiver  performance 
can  often  be  computed  directly  as  demonstrated  by  the  many  derivations  of  such 
probabilities  found  in  the  literature  (  see  [Ref.  2:  pp.  228-234)  and  [Ref.  3:  pp.  200-212] 
for  example  ).  However,  from  the  point  of  view  of  the  recipient  of  the  digital 
information,  it  is  the  Bit  Error  Ratio  (BER)  that  becomes  important  in  digital 
communication  applications.  While  results  on  the  BER  performance  of  M-ary 
modulation  receivers  is  far  more  limited,  closed  form  expressions  for  the  BER  of  M- 
PSK  receivers  have  been  obtained  by  Lee  [Ref.  4],  with  similar  results  independently- 
derived  by  Tan  [Ref.  5|. 

This  thesis  is  devoted  to  analyzing  the  application  of  direct  bit  detection  methods 
which  allow  the  recovery  of  bits  in  the  received  symbol  individually,  in  differential 
PSK  digital  communication  systems.  For  analysis  purposes,  8-DPSK  and  16-DPSK 
systems  are  considered  and  their  performances  evaluated.  The  particular  question  to 
be  answered  is  whether  for  8-DPSK  and  16-DPSK  modulated  signals,  direct  bit 
detection  methods  provide  comparable  performance  to  conventional  phase  detection 
receivers,  especially  under  conditions  of  improper  receiver  carrier  phase 
synchronization. 

In  this  introductory  chapter,  background  information  on  M-PSK  communication 
systems  is  given.  Known  optimum  receiver  structures  are  presented,  and  Direct  Bit 
Detection  Receivers  (DBDR)  are  introduced.  In  Chapters  II  and  Ilia  presentation  of 
the  transmitter  and  receiver  structures  respectively  is  given  along  with  their  operational 
properties,  and  the  logic  implemented  in  the  DBDR  for  8-DPSK  and  16-DPSK 


modulated  signals  is  discussed  in  detail.  Chapter  IV  includes  the  noise  performance 
analyses  of  8-PSK  and  8-DPSK  receivers,  and  Chapter  V  presents  similar  receiver 
analyses  for  16-PSK  and  16-DPSK.  In  both  of  these  chapters,  receiver  local  oscillator 
phase  errors  have  been  incorporated  so  as  to  be  able  to  account  for  degrading  effects 
due  to  improper  carrier  phase  synchronization.  The  obtained  mathematical  expressions 
of  receiver  performance  are  evaluated  on  a  computer  and  graphically  presented  as 
plots  of  bit  error  ratio  as  a  function  of  signal  to  noise  ratio  (SNR)  and  phase  errors. 

A.  BACKGROUND 

Direct  bit  detection  methods  utilize  two  signal  processing  channels  in  the  receiver 
in  order  to  recover  the  digital  information  without  phase  angle  measurements. 
Moreover,  these  methods  exhibit  optimum  performance  in  a  BER  sense  when  carrier 
synchronization  is  completely  achieved.  Because  of  the  attractiveness  ana  Simplicity  of 
DBDRs,  some  studies  have  been  carried  out  in  this  area,  and  a  modified  tw'o-channel 
receiver  for  8-PSK  has  been  built  and  its  noise  performance  measured  by  Thompson 
[Ref.  6].  Furthermore,  noise  performance  evaluations  of  DBDR  for  8-PSK  have  been 
presented  by  Myers  and  Bukofzer  [Ref.  7],  while  similar  analyses  for  the  16-PSK  case 
can  be  found  in  Reference  8.  Advantages  of  direct  bit  detection  methods  are  well 
documented  in  these  referenced  works,  in  which  all  analyses  have  been  done  assuming 
no  receiver  local  oscillator  phase  error,  Gray  Code  bit  to  symbol  mapping  and  an 
additive  white  Gaussian  noise  interference  model.  We  present  here  8-PSK,  8-DPSK, 
16-PSK  and  16-DPSK  communication  systems  under  similar  assumptions,  except  now 
the  receiver  has  an  assumed  known  local  oscillator  phase  error. 

B.  THEORY  OF  M-PSK  COMMUNICATION  SYSTEMS 

1.  Representation  of  M-PSK  Signals 

For  M-PSK  modulation  with  equal  signal  energies,  a  convenient 
representation  [Refs.  3:  p.p.  192-205]  of  the  signal  set  is  given  by 

Sj(t)  =  V2E7Tscos[27tfot+27t(i-l)/M],  i- 1,2 . M  ,  0  St  *Ts  (1.1) 


Es=  j  S;2(t)dt  ,  i—  1,2,. .,8 


Ts  =  symbol  duration 
f0  =  carrier  frequency  (Hz.) 
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A  suitable  orthonormal  signal  set  for  the  representation  of  the  above  signals  is 


<pj(t)=  cos(2Jtfot) 

<p2(t)=  V2?TS  sin(27tfot)  ,  0  S  t  ^  Ts  (1.3) 

so  that  by  using  trigonometric  identities,  s;  (t)  can  be  expanded  in  terms  of  <p.(t)  and 
<p-,(t).  That  is, 

Sj(t)  =  <p4(t)  cos[  2it(i-l)/M  ]  -  <f>  2(t)  Sint  2ic(i-l)/M  ]  (1.4) 

where  i=  1,2 . ,M  . 

A  plot  of  the  signal  constellation  for  M-PSK  is  shown  in  Figure  1.1  for 
various  values  of  M,  where  the  coordinates  of  the  i^  signal  vector  representing  the 
signal  s;(t)  are 

{V2ES  cos[2rc(i-l)/M] }  ,-  {  V2ES  sin(27t(i-l)/M] }  (1.5) 

For  M-DPSK,  the  signal  representations  are  exactly  same  as  those  already  presented, 
except  that  the  digital  information  is  transmitted  in  terms  of  phase  differences  between 
consecutive  signal  transmissions,  rather  than  absolute  signal  phases. 

2.  Gray  Coding  for  M-PSK  Signals 

Before  proceeding  with  the  Gray  coding  for  M-PSK  signals,  we  introduce 
some  definitions  which  will  henceforth  be  used 

1.  Symbol  :  Transmitted  or  received  signal  in  sinusoidal  form  (s;(t)  where 
i=  1,2,. ..M)  which  is  related  to  the  symbol  state. 

2.  Symbol  state  :  The  binary  digit  assignment  to  the  symbols  which  always 
consists  of  a  block  of  log2M  bits  in  length)^  where  i=  1,2,.. .,M.) 

3.  Data  state  :  Block  of  log2M  bits  corresponding  to  the  data  that  must  be 

delivered  to  the  intended  user,  and  which  is  processed  at  the  encoder  input  to 
determine  the  symbol  to  be  transmitted.  (Qi  where  i=  1,2 . M.) 

The  relationship  between  symbol  error  probability  and  bit  error  probability  in 
the  corresponding  log2M  bit  groups  depends  upon  the  assignment  of  bits  to  the  symbol 
states.  The  preferred  bit  to  symbol  state  mapping  in  most  cases  is  the  Gray  Code,  in 
which  the  symbol  states  which  correspond  to  adjacent  phases  dilTer  in  only  one  bit 
position.  This  is  because  in  the  coherent  demodulation  process  the  most  likely  error 


involves  mistaking  a  correct  symbol  with  its  nearest  neighbor  and  so  such  symbol 
errors  translate  into  single  bit  errors. 

In  Direct  Bit  Detection  Receivers  (DBDR),  Gray  code  assignment  to  the 
symbol  states  are  selected  for  the  above  given  reason.  In  this  thesis,  a  modified  Gray 
code  is  used,  which  for  comparison  purposes  is  shown  along  with  the  conventional 


TABLE  1 

COMPARISON  OF  CONVENTIONAL  AND  MODIFIED  GRAY 

CODING  SCHEMES 


Svmbol 

States 

Grav 

Coding 

Modified  Grav 
Coding 

8-PSK 

16-PSK 

8-PSK 

16-PSK 

Si 

000 

0000 

111 

mi 

S2 

001 

0001 

101 

1101 

S3 

Oil 

0011 

001 

1001 

S4 

010 

0010 

Oil 

1011 

S5 

110 

0110 

010 

0011 

S6 

111 

0111 

000 

0001 

S7 

101 

0101 

100 

0101 

s8 

100 

0100 

110 

0111 

S9 

1100 

0110 

Sio 

1101 

0100 

fin 

1111 

0000 

fi12 

1110 

0010 

Sl3 

1010 

1010 

S14 

1011 

1000 

fil5 

1001 

1100 

fi16 

1000 

1110 

Grav  Code  in  Table  1. 


Using  the  modified  Gray  Code  on  Table  1,  the  vectors  representing  the 
symbols  with  components  along  the  (pj  and  (p2  axes  have  been  shown  (see  Figure  1.2) 
such  that  Tj  is  the  symbol  component  along  the  axis  and  r,  is  the  symbol 
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component  along  the  <p,  axis.  Therefore  for  M  =  8,  (see  Figure  1.2.a)  in  order  to 
recover  Most  Significant  Bit  (LB)  of  the  transmitted  symbols  correctly,  r1  >  0  whenever 
the  symbol  states  S[y  Sj  and  g8  are  transmitted  and  r(  <0  whenever  the  symbol 
states  g3,  g4,  g-  and  g6  are  transmitted.  The  Middle  Bit  (MB)  of  the  transmitted 
symbols  can  be  recovered  correctly  if  (rj)2  -  (r,)2  >  0  whenever  the  symbol  states  gj, 
g4,  g5  and  g8  are  transmitted  and  (rj)2  -  (r2)2  <  0  whenever  the  symbol  states  g2,  £3. 
g6  and  g7  are  transmitted.  The  Least  Significant  Bit  (RB)  is  recovered  correctly  if 
r,>0  whenever  the  symbol  states  gj,  g2,  g3,  and  g4  are  transmitted  and  r, <0 
whenever  the  symbol  states  g5,  g6,  g7,  and  g8  are  transmitted. 

Similar  logic  has  been  set  up  for  the  M=  16  case  (see  Figure  1.2.b)  but  now 
each  symbol  is  related  to  4  bits  in  a  symbol  state.  Therefore  for  this  case,  4  bits  must 
be  recovered  from  each  received  symbol,  where  the  bits  in  a  symbol  state  are  labeled 
Most  Significant  Bit  (LB),  Middle  Bit  Right  (MBR),  Middle  Bit  Left  (MBL),  and  Least 
Significant  Bit  (RB).  For  example  from  Figure  1.2.b,  g7  is  0101,  corresponding  to 
LB  =  0,  MBL=1,  MBR  =  0  and  RB=1.  In  order  to  recover  the  LB  correctly,  it  is 
necessary  that  ^>0  for  transmitted  symbol  states  gp  g2,  g3,  g4,  g13,  g14,  g15,  and 
g16  and  that  rt  <  0  for  transmitted  symbol  states  g5,  g6,  g7,  g8,  g9,  g1Q,  gn,  and  g,2. 
If  r,>0  whenever  the  symbol  states  gj,  g2,  g3,  g4,  g5,  g6,  g7,  and  gg  are  transmitted 
and  r2<0  whenever  the  symbol  states  g9,  g10,  gH,  g12,  g13,  g14,  g15,  and  g16  are 
transmitted,  the  RB  will  be  recovered  correctly.  If  (rj)2>(r2)2  for  transmission  of 
symbol  states  gt,  g2,  g7,  g8,  g9,  g1Q,  g15,  and  g16,  and  (r3)2 < (r2)2  for  transmission  of 
symbol  states  g3,  g4,  g5,  g6,  g31,  g12,  g13,  and  g14,  the  MBL  will  be  recovered 
correctly.  Finally,  if  [  (rj)2  -  (r2)2  ]2  -  (2  r3 r2)2  >  0  upon  transmission  of  the  symbol 
states  gp  g4  ,  g5,  g8,  g9,  g12,  gn,  and  g16,  and  [  (rj)2  -  (r2)2  ]2  -  (2  r^)2  <  0  upon 
transmission  of  the  symbol  states  on  the  states  g„  g3,  g6,  g7,  g10,  gu,  g14,  and  g15, 
then  the  MBR  will  be  recovered  correctly. 

In  subsequent  chapters,  the  receiver  logic  and  its  performance  will  be 
explained  in  detail.  Additionally,  application  of  this  methodology  will  be  presented  in 
so  far  as  differential  PSK  modulation  is  concerned.  Before  doing  so  however,  a 
summary  of  known  results  on  the  performance  of  conventional  and  DBDR's  for  M- 
PSK  and  M-DPSK  is  presented  next. 

3.  Known  Results  on  the  Performance  of  M-PSK  and  M-DPSK  Receivers 

In  this  section,  known  results  on  the  receiver  performance  for  M-PSK  and 
DPSK  modulated  signals  are  presented  for  two  cases,  namely  conventional  receivers 
and  Direct  Bit  Detection  Receivers  (DBDR). 


a.  Conventional  M-PSK  and  M-DPSK  Receivers 

Coherent  demodulation  of  M-PSK  implies  in  principle  the  use  of  M  signal 
processing  channels  in  the  receiver,  as  shown  in  Figure  1.3  (from  Ha  [Ref.  10)  )  where 
different  forms  of  the  optimum  receiver  structure  of  M-PSK  signaling  have  been  given. 
Performance  analysis  results  for  optimum  receivers  are  well-documented  in  terms  of 
symbol  error  rate  as  a  function  of  signal  to  noise  ratio  [Refs.  2,9],  while  more  recent 
analyses  have  focused  on  the  receiver  probability  of  bit  error  (see  [Ref.  4]  and 
[Ref.  3:  p.p.  204-207  ]  for  example).  Symbol  error  probabilities  can  be  computed  from 
knowledge  of  the  channel  and  signal  characteristics.  However  for  transmission  of 
binary  data,  when  analyzing  system  performances  as  a  function  of  different  levels  of 
modulation,  the  bit  error  probability  rather  than  the  symbol  error  probability  is  of 
interest,  as  previous';,  .cplained.  In  Prabhu  [Ref.  1:  p.  198],  the  bit  error  probability 
has  been  presented  only  for  the  case  of  M-ary  orthogonal  signal  sets.  However,  since 
the  M-PSK  signal  set  is  not  orthogonal,  those  results  are  not  applicable  here.  The 
closed  form  expression  for  the  BER  of  M-PSK  has  been  determined,  when  a  Gray 
Code  bit  mapping  is  used  [Ref.2:  p.  198).  The  results  can  be  summarized  as  follows 

n 

Pb(M)  =  I(k/n)Pk  (1.6) 

k  =  1 

where 

M  =  Number  of  distinct  signal  waveform, 

Pk  =  Probability  of  k  bits  in  error  in  a  received  n  bit  data 
block  [Ref.  5:  pp.24-32], 
n  =  log2  M  , 

Pb(M)  =  Bit  error  probability. 

A  derivation  of  Equation  1.6  and  its  numerical  evaluation  (see  Figure  1.4  )  has  been 
carried  out  by  Tan  [Ref.  5:  p. 3 3].  For  M-PSK  signalling,  receiver  performance  analyses 
in  terms  of  bit  error  rate  have  also  been  presented  by  Lee  [Ref.  4:  p.49 1 J,  and  plotted  in 
Figure  1.5. 

b.  Direct  Bit  Detection  Receivers  ( DBDR's ) 

DBDR  structures  for  M-PSK  have  been  derived  for  various  values  of  M 
[Refs.  6,7,8).  As  previously  indicated,  this  method  is  particularly  well  suited  for  M-PSK 
communication  systems.  The  DBDR  structures  for  8-PSK  (see  Figure  1.6)  and  16-PSK 
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(see  Figure  1.8)  have  been  presented  and  their  performances  have  been  derived  in 
References  6,7,S.  The  BER  is  quantified  by  the  probabilities  of  a  bit  error,  denoted  PE. 

For  S-PSK,  PE  has  been  derived  by  Myers  and  Bukofzer  [Ref.  7:  p.  13] 

where 

PE=  1,3  {  Q[  2  VSNR  cosa  ]  +  Q[  2  VSNR  sina  ]  + 

+  Q[  2  v'SNR  (  cosa  +  sina)  ]  +  Q[  2  V^NR  (  cosa-  sina)  ]  - 
-  2  Q[  2  VSNR  (  cosa  -  sina)  ]  .  Q[  2  v'SNR  (  cosa+  sina)  ]  }  (1.7) 

where 

SNR  =  E,:2N0-3Eb.*2N0 
Es  =  Symbol  energy.' 

Ng=  Noise  power  spectral  density  lc.el  (one  sided) 
a  =  Basic  signal  phase  angle  (see  Figure  1.2). 

and 

oo 

Q(x)  =  J  [  1/V2*]  exp(-u2;2)  du  (1.8) 

x 

The  derivation  of  Equation  1.7  assumes  an  AWGN  channel,  Gray  coded  symbols  and 
equal  a  priori  bit  probabilities  (i.e.  ,  1/2).  Receiver  performance  in  BER  sense  is 
plotted  in  Figure  1.7  for  various  values  of  a  as  a  function  of  SNR. 

For  the  16-PSK  case,  probability  of  bit  error  is  derived  by  Bukofzer 
[Ref.  8]  and  the  resulting  receiver  performance  is  plotted  in  Figure  1.9  for  various 
values  of  the  angles  a  and  P,  where  P  is  shown  on  Figure  1.2.b. 

In  both  analyses,  it  has  been  assumed  that  the  receiver's  local  oscillator 
does  not  have  any  phase  error.  In  this  research,  similar  analyses  to  those  conducted 
in  References  7  and  8  have  been  carried  out  by  considering  the  effect  of  a  receiver 
local  oscillator  phase  error  for  8-PSK  and  16-PSK  cases,  as  well  as  the  noise 
performance  of  8-DPSK  and  16-DPSK  DBDR,  presented  as  BER  as  a  function  of 
SNR. 
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Figure  1.7  PE  Versus  SNR  (dB)  for  S-PSK  for  Various  Values  of  a. 


Figure  1.8  Block  Diagram  for  16-PSK.  DBDR. 
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II.  TRANSMITTER 


The  generation  of  differential  PSK  signals  can  be  performed  by  pre-processing 
the  data  bits.  This  chapter  is  devoted  to  analyzing  the  M-DPSK  transmitter  structure. 
One  general  form  of  the  transmitter  in  block  diagram  form  for  DPSK  communication 
systems  is  shown  in  Figure  2.1.  Explanation  of  the  general  elements  in  Figure  2.1  will 
be  included  along  with  the  derivation  of  the  transmitter  structure  for  8-DPSK  and 
1 6-DPS K  modulated  signals. 


Figure  2.1  DPSK  Transmitter  Block  Diagram. 


A.  TRANSMITTER  STRUCTURE  FOR  8-DPSK  AND  16-DPSK  MODULATED 

SIGNALS 

The  main  task  in  generating  the  DPSK  transmitter  involves  the  design  of  the 
encoder,  which  depends  on  the  algorithm  used  to  relate  data  states  to  the  symbol  states 
prior  to  transmission  as  a  sinusoidal  waveform.  Basically,  the  transmitted  symbol  is 
the  quantized  equivalent  of  the  angle  between  the  successively  transmitted  symbols. 
For  reasons  previously  explained,  a  Gray  coding  scheme  is  utilized  in  the  symbol  state 
to  symbol  mapping  for  the  two  special  cases  of  M-DPSK  signalling  considered  here, 
namely  M  *  8  and  M  =  16. 


The  symbol  differential  phase  relative  to  the  data  state  to  be  transmitted  is  given 
on  Table  2  for  8-DPSK  and  on  Table  3  for  16-DPSK  modulated  signals.  In  both 
Tables,  it  is  assumed  that  an  equiphase  arrangement  for  the  transmitted  symbols  is 


As  mentioned  earlier,  the  main  purpose  of  the  encoder  is  to  generate  a  mapping 
algorithm  whereby  the  phase  difference  of  consecutive  symbol  outputs  is  related  to  the 
input  data  in  accordance  with  the  information  in  Table  2  and  Table  3.  In  order  to 
achieve  successive  signal  transmissions  the  following  algorithm  is  implemented  in  the 
transmitter 

1.  Store  the  present  data  state  (£J;,  i=  1,2..., M)  and  the  previous  symbol  state 
(§.,i=l,2,3..,M). 

2.  From  the  Table  2  or  Table  3,  find  the  corresponding  differential  phase  value 
associated  with  the  present  data  state  as  Na,  where  N  is  an  integer  with  values 
of()<N<M-l. 

3.  Generate  a  symbol  state,  such  that  the  corresponding  symbol  will  have  a  phase 
angle  Na  degrees  larger  (  in  additive  sense,  modulo  2 n)  than  the  previously 
transmitted  svmbol. 
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For  example,  assuming  the  previous  symbol  state  is  S7  (0101)  and  present  data  state  is 
Dg  (0110)  which  corresponds  to  the  differential  phase  value  of  4a  (see  Table  3.)-  From 
Figure  1.2.b,  proceed  4a  degrres  counterclockwise  from  £7,  resulting  in  the  symbol 
state  is  g14  (1111)  will  be  the  symbol  output  from  the  encoder. 

This  algorithm  can  be  shown  in  matrix  structure,  in  which  the  row  entries  and 
the  column  entries  of  the  matrix  are  the  symbol  states  and  data  states  respectively,  and 
the  elements  of  the  matrix  are  the  symbol  state  which  correspond  to  the  appropriate 
phase  difference  in  the  transmitted  symbol  (see  Figure  2.3  and  Figure  2.2).  For 


example,  under  the  same  assumptions  given  in  the  previous  paragraph,  by  using  Pg 
(0110)  as  a  row  entry  and  £,  as  a  column  entry  in  the  matrix  in  Figure  2.3,  present 
symbol  output  from  encoder  is  found  as  (1111)  at  the  intersection  of  the  row  and 
column  of  the  matrix. 
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Figure  2.2  Encoder  Matrix  for  8-DPSK  Signals. 


The  hardware  implementation  of  the  encoder  structure  can  be  accomplished  in 
several  forms.  A  practical  approach  to  building  the  encoder  can  be  achieved  by  the  use 
of  EPROMs  (Erasable  Programable  Read  Only  Memory).  Once  the  truthe  table  the 
programming  of  the  EPROM  is  a  trivial  task  to  accomplish.  Demultiplexers  can  also 
be  used  to  build  the  encoder  circuitry:  however,  because  such  a  device  would  be  hard 
wired,  it  would  not  have  the  reprogrammability  of  the  EPROM  implementation. 

In  this  research,  the  8-DPSK  and  16-PSK  communication  systems  are  evaluated 
using  a  computer.  Appendix  A  includes  the  source  file  which  is  used  to  generate  the 
possible  outputs  of  an  M-PSK  and  an  M-DPSK  communication  system. 

The  carrier  modulator  uses  the  encoder  outputs  (symbol  states)  as  an  input,  and 
generates  the  symbols  (signals)  suitable  for  transmission.  The  signal  transmission 
involves  sinusoidal  voltages  whose  phase  is  related  to  the  symbol  states  (encoder 
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outputs)  and  of  carrier  frequency  suitable  for  the  transmission  channel.  The  signal 
generation  can  be  achieved  via  a  multiplexer,  whose  structure  is  shown  in  Figure  2.4. 
For  this  kind  of  implementation.  (M-l)  phase  shift  elements  are  needed  with  set  phase 
shift  values.  Another  approach  would  involve  an  implementation  using  a  D  A 
converter  along  with  a  phase  modulator  which  can  be  obtained  with  a  VCO  (see 
[Ref.  11:  pp.269-27 1J. 


Figure  2.4  Carner  Modulator  Structure  for  M-DPSK  Signaling. 


B.  COMPUTER  GENERATION  OF  THE  TRANSMITITR  OUTPUTS 

For  computer  evaluation  of  the  M-ary  DPSK  communication  system  ,  a 
FORTRAN  program  called  BELIZ  has  been  written  (see  Appendix  A).  The  program 
generates  all  possible  data  state  and  symbol  state  combinations  for  M-DPSK.  and 
produces  the  corresponding  symbol  states  outputs.  The  program  outputs  are  given  in 
Appendix  B  in  tabular  form  along  with  the  transmitter  matrix  of  Figure  2.3  and  Figure 
2.2,  specified  in  decimal  form. 


III.  RECEIVER 


A.  RECEIVER  STRUCTURE  FOR  S-DPSK  AND  16-DPSK  MODULATED 

SIGNALS 

The  introductory  chapter  presents  an  overview  of  the  structure  and  performance 
of  optimum  receivers  for  M-ary  signalling.  The  focus  of  this  chapter  is  on  DBDR's. 
As  mentioned  in  Chapter  II,  the  symbol  state  generation  based  on  the  data  state  and 
the  previous  symbol  state  is  the  most  important  task,  performed  by  the  encoder. 
Therefore  in  the  receiver’s  decoder,  a  similar  algorithm  must  be  performed  in  reversed 
sens*?,  that  is,  from  two  successive  symbol  states,  the  corresponding  data  state  must  be 
estimated.  A  generic  form  of  the  receiver  in  block  diagram  structure  is  shown  in  Figure 
3.1. 


Figure  3.1  General  receiver  block  diagram  for  DPSK  signals. 

The  algorithm  for  recovering  the  transmitted  data  is  as  follows 

1.  Recover  the  symbol  state  from  the  received  symbol  using  direct  bit  detection 
techniques. 

2.  Store  the  past  symbol  state  Shi  =  1,2.. ..M. 

3.  Compare  the  present  symbol  state  with  the  past  symbol  state  to  determine  the 
incremental  "  phase  angle  "  between  the  two. 

4.  From  Tabic  2  or  3  find  the  corresponding  data  state  that  gave  rise 

to  the  "phase  angle"  between  the  two  consecutive  symbol  states. 


'w5 


For  example,  assuming  the  past  symbol  state  is  (£6)  0001  and  the  present  symbol  state 
is  £3  (1001)  from  Figure  1.2.b,  the  angle  between  and  £$  is  3a.  From  Table  3.,  the 
corresponding  Gray  Code  representation  for  3a  is  found  to  be  (0010)  ,  therefore  the 
recovered  data  state  will  be  Z)p- 

The  decoder  implementation  can  be  depicted  as  a  matrix  structure  in  a  manner 
similar  to  that  shown  for  the  transmitter  encoder  in  Chapter  II.  For  8-DPSK  and 
16-DPSK  modulated  signals,  the  decoder  matrices  are  given  in  Figure  3.2  and  Figure 
3.3  respectively.  For  example,  assuming  the  same  conditions  as  in  the  above 
paragraph,  by  using  S6  and  S3  as  a  column  and  row  entries  for  the  matrix  given  in 
Figure  3.3,  the  recovered  data  state  from  the  decoder  is  determined  as  (0010). 


SYMBOL 
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000 

Fiaurc  3.2  Decoder  matrix  for  S-DPSK  sienals. 


B.  COMPUTER  GENERATION  OF  RECEIVER  OUTPUTS 

Appendix  C  includes  the  results  of  the  computer  evaluation  of  the  M-ary  DPSK 
receiver.  In  addition  to  generating  the  decoder  matrix  in  decimal  form,  the  output  for 
this  program  includes  all  possible  combinations  of  the  decoder  outputs. 
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Figure  3.3  Decoder  matrix  for  16-DPSK  signais. 
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C.  RECEIVER  LOCAL  OSCILLATOR  PHASE  ERRORS 

The  introductory  chapter  summarizes  the  studies  done  on  M-PSK  DBDR 
structures,  assuming  that  the  receiver  local  oscillator  has  no  phase  error.  In  this 
research,  a  fixed  receiver  phase  error  has  been  assumed  in  the  performance  analyses 
conducted,  and  its  effect  on  receiver  BER  is  presented  in  Chapter  IV  and  Chapter  V. 
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IV.  PERFORMANCE  ANALYSES  OF  8-PSK  AND  8-DPSK  DBDR'S  WITH 
LOCAL  OSCILLATOR  PHASE  ERROR 

In  previous  chapters,  DBDR's  are  introduced  and  their  operational  properties  are 
explained.  In  this  chapter,  performance  analyses  for  8-PSK  and  8-DPSK  will  be 
carried  out  in  detail. 

A.  8-PSK  DBDR  PERFORMANCE  ANALYSIS  IN  THE  PRESENCE  OF 
RECEIVER  LOCAL  OSCILLATOR  PHASE  ERROR 

The  possible  transmitted  signals  which  are  represented  as  vectors  in  Figure  1.2. a  , 
can  be  mathematically  expressed  as 

Sj(t)  =  V2E/T  sin  [  27tfQt  +  0.(t)  ]  ,  i-  1,2, 3..., 8  ,  0<  t  <T$  (4.1) 

where 


Es=  S  S;2(t)dt  ,  i=  1,2,. .,8  (4.2) 

0 

T  =  svmbol  duration 

s 

fQ  =  carrier  frequency  (Hz.) 

0.(t)  indicates  the  phase  of  the  transmitted  PSK  signal,  which  for  the  8-PSK  case 
becomes 


act)  - 


a  +  (i-1)  it,  4 
irt/4  -  a 


i=  1,3, 5, 7 
i=  2, 4, 6, 8 


0  <  t  <  Ts 


(4.3) 


The  angle  a  is  variable  as  shown  in  Figure  1.2. a;  however,  the  constraint 
0  <  a  <  n/4  must  be  observed.  Equation  4.1  can  be  written  by  using  a  suitable 
orthonormal  set  of  functions  as 


Sj(t)=  <pt(t)  7ES  cos0.(t)  +  <p,(t)  V/Es  sin0;(t) 


(4.4) 


i  —  1,2 S  ;  0<t<T$ 

where  cpj(t)  and  q>-,(t)  are  the  orthonormal  function  set  which  by  applying  the 

Gram-Schmitt  orthonormalization  procedure 

[Ref.  3:  p.l9S] 
can  be  shown  to  be 

cpj(t)  =  V2  T  sin  2rtfgt 

tp2(t)  =  VFTcos  2nf0t  ,  0<  t  <TS  (4.5) 

The  received  signal  depends  upon  the  channel  model  and  the  transmitted  signal 
as  specified  by  Equation  4.1.  In  this  research,  an  Additive  White  Gaussian  Noise 
(AWGX)  interference  model  has  been  assumed  so  that  the  received  signal  can  be 
written  as 

r(t)=S;(t)  +  n(t)  (4.6) 

where  n(t)  is  a  sample  function  of  a  white  Gaussian  noise  process  having  a  constant 
(one  sided)  power  spectral  density  function  equal  to  NQ  watts; Hz.  The  receiver 
processes  r(t)  over  an  interval  of  duration  T$  seconds  and  a  decision  is  made  at  the  end 
of  the  interval.  This  procedure  continues  in  a  repetitive  fashion  as  time  progresses. 

Assuming  all  signals  are  equally  likely  to  be  transmitted,  the  probability  of  a  bit 
being  correct  for  M-PSK  signal  transmission  is  written  as 

M 

Pr{  bit  correct  }  =  1,'M  £  Pr{  bit  correct  !  s;(t)  transmitted  }  (4.7) 

i=  1 


mr- 


tv 


K 


V  vv  V  V 


The  calculations  for  the  conditional  probabilities  (Pr{  h:t  correct  '  Sj(t)  transmitted  }, 

i=  1,2 . M  )  have  been  determined  by  Myers  and  Bukofzer  [Ref.  7:  pp.  6-12],  for  the 

8-PSK  case  with  no  receiver  local  oscillator  phase  error.  In  the  sequel,  the  same 
analysis  methods  are  utilized  to  obtain  the  receiver  BER  in  the  presence  of  a  local 
oscillator  phase  error.  A  local  oscillator  phase  error  is  modeled  by  the  angle  c  in  the 
unmodulated  receiver  carriers  <pj(t)  and  <p,(t)  of  Figure  1.6.  These  are  denoted  <p t '( t ) 
and  <p2'(t),  namely 
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<pj'(t)  =  J2,  T  sin  (  27tfQt  +  z  ) 

(p;'(t)  =  VTT  cos  (  27Tf0t  +  e  ),  0  <  t  <  TS  (4.S) 

Assuming  signal  s;(t)  was  transmitted,  this  produces  outputs  rk  .  from  the  receiver's 
correlators,  where 


rk,i  =  J  t  +  n(c)  1  %  dt  ’  i=  1*2»-8  *  k.=  1,2  (4.9 

0 

For  example  when  Sj(t)  is  transmitted,  from  Equation  4.4,  0.(t)  becomes  just  a  and 


rU  =  VHS  cos  (a-s)  +  nj 
r2  j  =  VEs  sin  (a-£)  +  n2 


(4.10) 


where 


nk  -  j"  n(t)  <pk'(t)  dt  .  k=  1,2 
0 

Observing  that 


E  [n,  n2]  =  E  [  J  n(t)  (pj'(t)  dt  x  j  n(t)  <p2'(t)  dt  ]=  0 
0  0 


(4.11) 


(4.12) 


This  means  that  nj  and  n,  are  uncorrelated  random  variables.  Since  n;  and  n,  are 
Gaussian  random  variables  (having  identical  variances)  [Ref.  7:  p.7],  they  are  also 
statistically  independent. 

The  goal  of  the  receiver  of  Figure  1.6  is  to  correctly  recover  the  LB,  the  MB.  and 
the  RB,  transmitted  via  the  signal  s;(t).  By  introducing  the  more  compact  notation 

Pr{  bit  correct  /  s;(t)  }  =  Pr(  bit  correct  S;(t)  transmitted  }  ,  i=  1,2. ..S  (4.13) 


CO 


the  probabilities  of  such  correct  bit  recovery  are  given  in  general  by 


Pr(  LB  correct  /  s;(t) }  = 


Pr(  r1(i>0  /  Sj(t) }  i=  1,2,8, 7 
Pr(  ri  ;  <  0  /  s;(t)  }  i=  3, 4, 5, 6 


(4.14) 


Pr f  rl  /  '  r2,i2  >  0  ''  si( 1  =  1  >4’5’ 


Pr{  MB  correct  /  Sj(t) }  = 


1,4, 5, 8 


(4.15) 


Pr{  rl,i2  ‘  r2,i2  <  0  '  Si(t)  J  i=  2’3’6>7 
Pr(  r,  .  >  0  /  s(t) }  i=  1,2, 3, 4 

— J1  i 

Pr{  r,  .  <0  /  s.(t)  )  i=  5, 6, 7, 8 

I 


Pr(  RB  correct  /  s^t)  }  = 


(4.16) 


From  Equations  4.8-4.16,  we  can  write 

Pr(  LB  correct  /  s^t) }  =  Pr(  rt  j  >  0  /  Sj(t) } 

=  Pr{  nj  >  -  cos(a-e)  /  Sj(t) } 
=  1  -  Q  [  J2  Es;N0  cos(a-e)  ] 

where 

oo 

Q(x)  =  J  [  1/-/2K  ]  exp(-u2/2)  du 


(4-17) 


(4.18) 


Similarly 


Pr(  RB  correct  /  Sj(t)  }  = 


8  Pr{  r2)1  >  0  /  Sj(t) } 

=  Pr{  n2>  -  sin^a-e) /  st(t) } 
=  1  -  Q  [  vTE7^0  sin(a-e)  ] 


(4.19) 


In  order  to  evaluate  the  conditional  probability  of  the  MB  being  correct,  the  following 
steps  are  required.  From  Equation  4.15 


Pr{  MB  correct  /  Sj( t)  }  = 


ri,i2-r2,i2  >0'si(t)) 

Pr«rl,i-  r2,i)-(rl,i+  r2,i>  >0/  S^)} 


(4.20) 


b 


f 

i 

t 
i 

\ 

By  defining  variables  yl  i  and  y2  ;  as 

I 

* 


then  from  Equation  4.9,  when  s1(t)  is  transmitted, 

Vj  j  =  V  Es  [  cos  (a-e)  +  sin  (a -£)  ]+  nt  +  n2 
y2  I  =  V  Es  [  cos  (a-e)  -  sin  (a-£)  ]+  nt  -  n2 
Clearly  v,  2  and  y,  .  are  Gaussian  random  variables  with  statistical  means  given  by 


E(  >'i  i)  =  VE^  [  cos  (a-£)  +  sin  (a-£)  ]  (4.24) 

E  t  >*2  i }  =  [  cos  (a-£)  -  sin  (a-£)  ]  (4.25) 

and  variances  given  by 

Var{yu}  =  Var{  y2>1}  =N0  (4.26) 

By  observing  that 

E  [  (ylfl-  E  [  ylfl])  (y2jl-  E  [  y2jl]) }  =  E  {  (nx  +  n2)  (nr  n2) } 

=  E  { (n,)2-(n2)2  }  =  0  (4.27) 


where  nt  and  n2  are  defined  in  Equation  4.11.  Since  yt  l  and  y2  j  are  uncorrelated 
Gaussian  random  variables,  they  are  also  statistically  independent.  Therefore, 
Equation  4.20  becomes 

Pr{  MB  correct  s1(t)  }  =  Pr{  Vj  j  y2  j  >  0  /  Sj(t)  ) 

=  Pr(  [  y’j  ! >  0  /  s,(t)  ]  and  [  y2>1  >  0  /  s,(t)  ]  }  + 

+  l  Yi,i  <0  '  si<l)  ^  and  l  >2,1  <0  '  si(t)  1  >  + 

=  Pr{  y1(1  >  0  /  Sj(t)  }  Pr(  y21  >  0  ,  Sj(t)  }  + 

+  Pr(  yM  <0  st(t)  }  Pr{  y2  l  <0  /  Sj(t)  } 


(4.21) 

(4.22) 

(4.23) 


(4.28) 


From  Equations  4.22  -  4.26,  we  have 


(4.29) 

(4.30) 


} 

I 
I 

Pr {  yu  >0  /  Sj(t) }  =  Q  (  E{Vj  j}  /  Vn0  ] 
pr(  y2  l  >0  /  s,(t) }  =  Q  [  E{y21}  /  3 

so  that  Equation  4.28  becomes 

Pr{  MB  correct  /  Sj(t) }  =  1  -  Q{  V  ES/NQ  [  cos(a-s)  4-  sin(a-£]  }  - 

-  Q{  V  Es/N0  [  cos(a-e)  -  sin(a-£] } 

-  2  Q{  J  Es,  N0  [  cos(a-e)  -  sin(a-e] }  . 

.  Q{  V  Es/N0  [  cos(a-s)  +  sin(a-e] }  (4.31) 

Assuming  now  transmission  of  s2(t),  a  similar  analysis  procedure  is  performed 
but  now  Equations  4.14  -  4.16  become 

02(t)  =  tr/2  -  a  (4.32) 

ri  2  =  sin  (a  +  e)  +  nj 

r2,2  =  VEjcos(a  +  £)  +  n2  (4.33) 

where  nt  and  n2  are  defined  by  Equation  4. 1 1.  From  Equations  4.14  and  4.33 

Pr(  LB  correct  /  s2(t) }  =  Pr{  Tj  2>  0  /  s2(t) } 

=  Pr{  n2>  -  sin(a  +  £)  /  s2(t) } 

=  1  -  Q  [  V 2  Es/N0  sin(a  +  e)  ]  (4.34) 

and 

Pr(  RB  correct  /  s2(t) }  =  Pr(  r2  2>  0  /  s2(t) } 

■  Pr{  n2>  -  vTT  cos(a  +  £)  /  s2(t) } 

=  1  -Q[  J2  Es/N0  cos(a  +  £)  ]  (4.35) 


For  calculations  involving  the  MB,  one  additional  random  variable  needs  to  be  defined, 


Assuming  transmission  of  s,(t),  Equation  4.36  becomes 


f2,2  ‘  fl,2 


=  V  Es  cos  (a  +  £)  +  n,  -  V  Es  sin  (a  +  £)-nj 


(4.37) 


where 


E  {y3  2}  =  cos  (a -he)  -  y~E^  sin  (a  +  E)  =  E  {y,^} 


(4.38) 


Since 


Pr(  MB  correct  '  s2(t)  }  =  Pr{  (r1;1)2  -  (r2)2)2  <  0  /  s2(t) } 

=  Pft  <rl,l  *  r2,2>  <rI,l  +  r2,2)<0  '  S2W  > 

bv  using  Equations  4.21  and  4.36 

Pr(  MB  correct /  s2(t) }  =  Pr(  y,j2  y3  2  <  0  /  s2(t) } 

Furthermore 

E  {  (y1)2-  E  (  >’1,2))  (y^X  E  {  y3,2J)  J  =  E  Uni  +  n2)  (n,-  n2) } 

=  E  {  (nj)2-(n2)2  }  =  0 


(4.39) 


(4.40) 


(4.41) 


so  that  yl  2  and  y3  2  are  statistically  independent  Gaussian  random  variables,  and 
therefore  Equation  4.40  becomes 


Pr(  MB  correct  /  s2(t) }  =  Pr{  y3>2  y1>2  <  0  /  s2(t)  } 

•  Pr(  lyi,2>0  /  M1)!  and[y3  2>0,'  s2(t)  ]  }  + 

+  Pr{  [  y1)2<0  '  s2(l)  1  and  [  ^3.2<0  '  s2^  ]  }  + 
=  Pr{  yli2>0  !  S2(t) )  Pr{  y3>2>0  s,(t)  }  + 

+  Pr{  yL  2<0  ,  s2(t)  }  Pr{  y3,<0  s2(t)  } 


(4.42) 


By  using  a  method  similar  to  that  used  in  the  derivation  of  Equations  4.29  and  4.30. 
for  Vj  2  and  y3  2  ,  Equation  4.42  becomes 


WWW? 


I 


£ 


K 


K 


£ 


•5 

'1 


It 


V. 


Pr(  MB  correct  s2(t) }  =  1  -  Q(  J  Es  N'0  [  cos(a  +  £)  +  sin(a  +  £] 


-  Q(  V  E$  Ng  [  cos(a  +  e)  -  sin(a  +  £]  } 


-  2  Q{  J  Es  Nq  [  cos(a  +  £)  +  sin(a  +  £]  } 


.  Q(  J  E  Ng  [  cos(a  +  £)  4-  sin(a  +  £] 


(4.43) 


Analyses  can  be  carried  out  for  probability  of  correct  reception  assuming  even7 
possible  signal  transmission.  This  is  however  not  necessary  because  of  inherent 
symmetries  in  the  problem  .  In  conditions  where  symmetries  are  not  present,  previous 
methodologies  can  be  used  to  obtain  results  as  follows.  Assuming  transmission  of 
s3(t),  Equations  4.3  and  4.9  become 


03(t)  =  n;  2 


a 


Pr{  LB  correct  /  s3(t) }  =  Pr(  rt  3  <  0  /  s3(t)  } 

=  Pr{  >  -  VES  sin  (a-£)  /  s3(t) } 
=  1  -  Q[  J2  Es/Ng  sin(a-£)  ] 


and  from  Equation  4.16 
Pr(  RB  correct  /  s3(t) } 


=  Pr{  r2  3  <  0  /  s3(t) } 

=  Pr{  n2>  -  VES  cos  (a-E)  /  s3(t) } 
=  1  -  Q[  Jl  Es/N0  cos(a  +  e)  ] 


Assuming  transmission  of  s,(t),  Equations  4.3  and  4.4  become 


04(t)  =  k-  a 

r j  4  =  7ES  cos  (7t-a-£)  +  nt 
=  -  vE^  cos  (a  +  £)  +  nj 


44 


(4.44) 


0 

rj  3  =  vEs  cos  (a-£  +  7r/2)  +  nj 

(4.45) 

ijj 

=  -  VES  sin  (a-£)  +  nt 

r2,3  =  sin  (a-£  +  jr/2)  +  n2 

(4.46) 

1 

=  -  v'Ej  cos  (a-£)  +  n2 

V 

i 

so  that  from  Equation  4.14 

(4.47) 


(4.48) 


(4.49) 

(4.50) 


r2  4  =  v Es  sin  (Ji-a-e)  +  rt, 
=  •  vlT  sin  (a  +  £)  +  n. 


so  that  from  Equation  4.14 


Pr{  LB  correct  !  s4(t)  }  = 


:  Pr(  ri,4  <  0  ■'  s4(0  ) 

=  Pr{  nj  >  -  sin  (cH-£) !  s4(t)  } 
=  1  -  Q[  *J2  Es  Nq  cos  (a  +  £)  ] 


and  from  Equation  4.16 


Pr(  RB  correct  /  s4(t)  }  = 


=  Pr{  r2>4  >  0  ,  s4(t)  } 

=  Pr{  n2>  -  jEs  sin  (a  +  £)  /  s4(t) } 
=  1  -Q[  VTE7N0sin(a  +  £)] 


Also  from  Equation  4.15 

Pr(  MB  correct  /  s4(t) }  =  Pr{  (rlj4)2  -  (r2j4)2>  0  /  s4(t) } 

=  Pr(  <rl,4  '  ‘V  <rl,4  +  r2,4>>0  ''S4<1>> 

Using  Equations  4.21  and  4.36  yields 


'1 ,4  =  rl,4  +  r2,4 
■2,4  =  rl,4  r2,4 


and  using  a  method  similar  to  that  used  in  the  derivation  of  Equation  4.27  results  in 

Pr(  MB  correct  /  s4(t) }  =  Pr{  ylf4  y2<4  >  0  /  s4(t) } 

=  Pr(  [  yli4:>0  /  s4(t)  ]  and  [  y2  4>0  /  s4(t)  ]  }  + 

+  Pr{  [  y14  <  0  /  s4(t)  ]  and  [  y2<4  <  0  ;  s4(t)  ]  }  + 

=  Pr{  y14>0  -  s4(t) }  Pr{  y2,4>0  s4(t)  }  + 

+  Pr{  ylj4<0  s4(t)  }  Pr{  y24<0  /  s4(t)  }  (4.5' 

which  becomes 


i*.  «*.  iC 


Pr{  MB  correct  /  s4(t) }  =  1  -  Q  [  v'2  E$  NQ  [  cos(a  +  £)  -  sin(a  +  £)  ]  }  . 

-  Q  {  E$  Nq  [  cos(a  4  e)  +  sin(a  4-  z)  ]}  + 

+  2  Q  (  v'2  Es  Nq  [  cos(a  +  £)  -  sin(a  +  £)  ]  }  . 

.  Q  {  v2  Es  N0  [  cos(a  +  £)  4-  sin(a  +  E)  ]}  (4.58) 


For  transmission  of  s5(t),  the  pertinent  equations  are 


05(t)  =  71+  a 

(4.59) 

rj  5  =  v^Es  cos  (7t  +  a-£)  +  nt 

(4.60) 

=  -  VES  cos  (a-£)  +  nj 

r2  5  =  sin  (7t  +  a-£)  +  n2 

(4.61) 

=  -  v"Es  sin  (a-£)  +  n2 

so  that 

Pr{  MB  correct  /  s5(t) }  =  Pr(  (r15)2  -  (r2  S)2>  0  /  s$(t)  }  (4.62) 

Pr(  (r1)5  -  r2  5)  (r15  4-  r2(5)>0  /  s5(t) } 

which  becomes 

Pr{  MB  correct  '  s$(t)  }  =  Pr{  Yi  5  5  >  0  ss(l)  ) 

=  Pr{  [  yii5>  0  /  s5(t)  ]  and  [  y2  5>  0  /  s5(t)  ]  }  + 

+  Pr(  t  >'i,5<0 ' 5 5(0 )  and  l  y2,5<0-'  s5(0  ]  )  + 

=  Pr(  y1?5>0  /  s5(t) }  Pr{  y2?5 >  0  /  s5(t)  }  + 

+  Pr{  yi’5<0  /  s5(t) }  Pr{  y2’5<0  /  s5(t)  }  (4.57) 

by  observing  that  Vj  5  and  y2  5  are  statistically  independent  random  variables  making 
it  possible  to  apply  the  similar  procedure  shown  above.  The  result  becomes 

Pr(  MB  correct /  s$(t)  }  =  1  -  Q  {  J2  ES/NQ  [  cos(<z-£)  +  sin(a-£)  ]  }  - 

-  Q  {  J2  Es/N0  [  cos(a-£)  -  sin(a-£)  ]}  + 

+  2  Q  {  y/2  Es/Ng  [  cos(a-£)  +  sin(a-£) )  }  . 

.  Q  {  J2  E  :N0  [  cos(a-E)  -  sin(a-E)  ]} 


(4.64) 


I  M  XI '  !.• 
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As  indicated  before,  because  of  the  inherent  symmetry  in  the  problem,  the 
following  equations  hold  for  conditional  LB  correct  probabilities  (see  Figure  1 . 2. b  ) 


Pr(  LB  correct  ,  Sj(t)  }  : 
Pr(  LB  correct  '  s,(t)  }  : 
Pr(  LB  correct  s3(t)  }  : 
Pr(  LB  correct  /  s4(t) }  : 

and  similarly  for  the  RB 

Pr(  RB  correct  '  Sj(t)  } 
Pr(  RB  correct  s,(t)  } 
Pr{  RB  correct  /  s,(t) } 
Pr[  RB  correct  /  s4(t)  } 

and  finally  for  the  MB 

Pr{  MB  correct  /  Sj(t)  } 
Pr(  MB  correct !  s2(t)  } 
Pr(  MB  correct  /  s4(t)  } 
Pr{  MB  correct  /  s5(t) } 


Pr(  LB  correct  !  s5(t)  } 

Pr(  LB  correct  s6(t)  } 

Pr(  LB  correct  s?(t)  } 

Prf  LB  correct  /  s8(t) } 


(4.65) 

(4.66) 

(4.67) 

(4.68) 


Pr{  RB  correct  ;  s5(t) } 
Pr {  RB  correct  /  s6(t)  } 
Pr(  RB  correct  /  s7(t)  } 
Pr(  RB  correct  /  sg(t)  } 


(4.69) 

(4.70) 

(4.71) 

(4.72) 


Pr{  MB  correct  /  s3(t) } 
Pr(  MB  correct  /  sg(t)  } 
Pr{  MB  correct  /  s6(t) } 
Pr(  MB  correct  /  s7(t)  } 


(4.75) 

(4.76) 

(4.77) 

(4.78) 


By  using  now  Equation  4.7  for  the  unconditional  probability  of  the  LB,  MB  and 
RB  being  correct,  the  above  results  yield 

Pr(  RB  correct  }  =  Pr{  LB  correct } 

=  1  -  1/4  [Q[  V/2E7N0cos(a-E)]  + 

+  Q[  %/2  Es/Nq  sin(a  +  £)  ]  + 

+  Q[  Es  Nq  cos(a  +  e)  ] 

+  Q[  vTE^Tq  sin(a-e)  ]  ]  (4.79) 


Pr(  MB  correct  }  =  1  -  1/2  [  Q  (  J2  Es,  NQ  [  cos(a-c)  +  sin(a-e)  ]  }  + 
+  Q  (  72  Es;N0  [  cos(a  +  e)  +  sin(a  +  £)]}  + 


.•  n-.nv.  .• 


W! 


+■  Q  !  v  -  E,  Ng  [  cosia-e)  -  sin(a-e)  ]  ]  + 

+  Q  :  v  -  Es  N(J  [  cosia  i-  o  -  sin(a-t-E)  ]}  ]  + 

~  Q  :  v  2  E.  Y,  [  cc -(<*-£)  -  sin(a-c)  ]  }  . 

Q  v  2  V.  Y,  |  cosia -ci  sinia-o  ]}  -i- 

*  Q  v  2  E  Nf)  |  cosia-ei  -  sinia-o  ]  )  . 


2  1  Y  '  com  a  -  ci  -  sinta  o 


(4.80) 


For  8-PSK.the  Signal  to  \oise  Ratio  'SNRi  has  been  defined  by  Equation  1.6,  so 
that  using  Equations  4  ~'i  -  4  mi.  -he  HI  R  denoted  PE  can  be  expressed  as  a  function 
of  SNR,  a  and  c  as  follows 


PE=  1/3  f  Pr(  LB  correct  -  pr  MB  correct  ,  ~  Pr‘  D  B  correct  ) 
=  1  -  1  6  {  Q|  2  v  SNR  com  a-E •  -  Oj  2  v  SN R  sin( a  -  c I  j  i- 
+•  Q[  2  v SNR  cosia  -  ci  ]  *  Q[  2  v  SNR  sinia-o  j 


comu-Ci  sima-ci 


+  Q  {  2  v;  SNR  [  cosi a  -  c i  *  smi a  *  £ )  ]  ) 
f  Q  {  2  v  SNR  [  cosia -ci  -  sinia-ci  j  ;  - 

+  Q  (  2  v  SNR  [  cosia  -  ci  -  sima  -  e )  ]  ;■ 

+  Q  {  2  v;  SNR  [  cosia-e)  -  sinia-o  ]  ). 

.  Q  (  2  J  SNR  [  cos(a  +  e)  +  sinta  +  £)  ] 

+  Q  {  2  -J  SNR  [  cos(a-c)  -  sinia-o  |  }. 

.  Q  {  2  V  SNR  [  cos(a  +  c)  -  sin(a  +  £)  j  } 


(4.81) 


From  Equation  4.81,  the  PE  can  be  computed  for  various  values  of  SNR,  a.  and 
£,  with  the  results  used  to  optimize  the  DBDR  performance.  For  £=  0°  (i.  e.  ,  no 

receiver  local  oscillator  phase  error),  values  of  PE  vs  SNR  have  been  plotted  (see 
Figure  4.1)  showing  complete  agreement  with  the  results  obtained  by  Myers  and 
Bukofzer  [Ref.  7:  p.  11]  Figure  4.1  includes  also  the  PE  for  various  values  of  £,  as  a 
function  of  SNR. 


8-DPSK  DBDR  PERFORMANCE  ANALYSIS  IN  THE  PRESENCE  OF 
RECEIVER  LOCAL  OSCILLATOR  PHASE  ERROR 


The  advantages  and  disadvantages  of  differential  encoding  and  decoding  have 
been  previously  explained  in  the  introductory  chapter.  In  this  section,  performance  of 
the  8-DPSK  DBDR  is  presented.  As  explained  earlier  the  8-DPSK  decoder  utilizes  the 
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Figure  4.1  PE  Versus  SNR  of  8-PSK  DBDR  for  Various  Values  of  e  (U 


LB,  MB  and  RB,  which  are  recovered  via  direct  bit  detection  methods,  and  performs  an 
inverse  mapping  from  which  the  original  data  stream  is  delivered  to  the  intended  user 
Before  proceeding  with  the  performance  calculations,  the  following  notation  will 
be  introduced.  The  probabilities  of  a  bit  being  correct  or  being  in  error  will  be  denoted 

Ck  =  Pr{  k  correct  }  ,  Ek  =  Pr{  k  error  }  (4.82) 

respectively,  where  k=  1  for  LB,  k=  2  for  MB,  k=  3  for  RB.  By  introducing  the 
notation  given  by  Equation  4.82,  the  analysis  for  the  performance  of  the  8-DPSK 
DBDR  can  be  expressed  in  more  compact  form.  For  example  E2  means  the 
probability  that  the  MB  is  in  error. 

In  order  to  obtain  the  probability  of  a  bit  being  correct  (with  differential 
decoding)  we  must  first  evaluate 

Pr{  RB  correct }  =  1/2  Pr{  RB  correct  /  RB=  1  }  + 

+  1/2  Pr{  RB  correct  /  RB=  0  }  (4.83) 

where  by  further  conditioning,  we  equivalently  have 

1  1 

Pr{  RB  correct/RB  =  m  }  =  £  £  Pr(  RB  correct  /  RB  =  m,MB=i,LB  =  j  }. 

i  =  0  j  =  0 

.  Pr{  MB  =  i,  LB  =  j  }  m  =  0.1  (4.84) 

Due  to  the  assumption  of  equally  likely  probabilities  of  having  0's  and  l's  in  the  data, 
Pr{  MB  =  i,  LB  =  ]  }  =  1/4  i,j=  0,1  (4.85) 

so  that  Equation  4.83  becomes 

1  1 

Pr{  RB  correct }  =  1,8  £  £  Pr{  RB  correct  /  RB  =  0.MB  =  i,LB  =  j  )  + 

i  =  0  j  =  0 


(4.S6) 


1  1 

4-  1/8  V  V  Pr(  RB  correct  /  RB=  1,MB  =  i,LB  =  j  } 
i  =  0  j  =  0 

The  DBDRs  recover  from  the  received  signal,  each  bit  in  the  data  state  (for  PSK 
modulation)  or  each  bit  in  the  symbol  state  (for  DPSK  modulation)  independently. 
This  property  of  the  DBDR  is  an  advantage  for  BER  analyses.  Because  certain 
consecutive  bit  errors  in  the  past  and  presently  recovered  symbol  states  do  not 
necessarily  result  in  errors  in  the  recovered  data  states,  these  bit  errors  in  successive 
symbol  states  must  be  considered.  That  is,  even  though  symbol  errors  may  occur,  the 
delivered  data  may  stillbe  correct.  Therefore,  these  type  errors  have  to  be  accounted 
for  as  resulting  in  a  correct  reception  and  considered  in  the  computation  of  the 
probability  of  a  bit  being  correct.  Figure  4.2  shows  one  example  of  the  situation  being 
discussed.  Here  the  data  state  (0,0,0)  is  assumed  to  occur  and  the  symbol  state  (0,0,0) 
is  assumed  to  have  been  previously  transmitted.  This  means  that  according  to  the 
encoder  matrix  shown  in  Chapter  II,  the  symbol  state  (0,0,0)  must  now  be  generated 
for  transmission.  Figure  4.2  further  shows  the  receptions  of  the  symbol  states  which 
even  though,  most  contain  errors,  produce  recovered  data  states  in  which  the  RB  is 
correctly  delivered.  Observe  that  there  are  32  such  combinations  which  yield  a  correct 
0  for  the  RB.  Therefore  notation  has  to  be  introduced  in  order  to  account  for  error 
patterns  of  the  above  type,  namely  Pr{  RB  correct,  (k,l,m)  /  RB  =  0,MB  =  i,LB  =  j  }, 
where  k.l,m=  0,1.  The  meaning  of  this  expression  is  the  probability  of  the  RB 
correctly  recovered  given  that  the  data  state  was  (j,i,0)  and  the  delivered  data  state  is 
(k,l,m)  (where  m=0  in  this  case,  otherwise  the  RB  could  not  be  correct  due  to  the 
transmission  of  the  RB  =  0).  Therefore 

Pr{  RB  correct  !  RB  =  0,MB  =  0,LB  =  0  }  = 

1  1 

=  V  £  Pr(  RB  correct. (k.1,0) ,  RB  =  0.MB  =  0,LB  =  0  }  (4.87) 

k=0  1=0 


Since  each  element  of  the  symbol  state  is  independently  recovered,  we  see  from  Figure 
4.2  that,  for  example 
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CORRECT  RECEPTION  COMBINATIONS  IN  PRESENCE  OF  ERROR 


INPUT  PREVIOUS  PRESENT  PREVIOUS  PRESENT 

TRANSMISSION  TRANSMISSION  RECEPTION  RECEPTION 


OUTPUT 


0  0  0  0  0  0 


0  0  0 


0  0  0 
0  0  I 

0  10 
Oil 
10  0 
1  0  1 
1  1  0 
111 
0  0  0 
0  0  1 
0  10 
0  1  1 
10  0 
1  0  1 
1  1  0 
1  1  1 
0  0  0 
0  0  1 
0  10 
0  1  1 
1  0  0 
10  1 
1  1  0 
1  1  1 
0  0  0 
0  0  1 
0  10 
Oil 
1  0  0 
10  1 
1  1  0 
1  1  1 


0  0  0 
0  0  1 
0  10 
Oil 
1  0  0 
10  1 
1  1  0 
1  1  1 
111 
0  0  0 
110 
10  0 
1  0  1 
0  1  0 
0  0  1 
0  1  1 
0  1  0 
1  0  1 
0  1  1 
0  0  1 
0  0  0 
1  1  1 
10  0 
1  1  0 
1  0  1 
1  0  0 
111 
1  1  0 
0  0  1 
0  0  0 
Oil 
0  10 


0  0  0 

0  0  0  i 

0  0  0' 

0  '  0  > 

0  0  0  ! 

0  0  0 

0  0  0 

0  0  0 

0  1  0  I 

0  10 
0  10 
0  10 
0  1  0 

0  1  0 

0  1  0 

0  10  1 
10  0; 
10  0 
1  0  o  I 

10  0 
10  0 
10  0 
10  0 
10  0 
1  1  0 

1  1  0 
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1  1  0 

110 
110 
1  1  0  i 

110  1 


Figure  4.2  Example  of  Error  Patterns  for  S-DPSK  System  Where  Data  State. 
Previous  Transmission,  and  Present  Transmission  are  000  with  RB  Correctlv  Recovered. 
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Pr{  RB  correct, (0,0,0)  /  RB  =  0,MB  =  0.LB  =  0  }  =  Cl2  C22  C3:  + 

+  Cl2  C22  E32 ■+•  Cl2  E22  C32+  Cl2  E22  E32+  El2  C22  C32  + 

+  El2  C22  E32+  El2  E22  C32+  El2  E22  E32  (4.38) 

From  Equations  4.79  and  4.82  -  4.86,  it  is  clear  that  Cl=  C3  and  EI=  E3. 
Therefore,  simplification  of  Equation  4.88  is  possible,  thus  yielding 

Pr{  RB  correct, (0,0,0)  RB  =  0,MB  =  0,LB  =  0  }  =  Cl4  C22+  2  Cl2  C22  El2  + 

+  Cl4  E22+  2  Cl2  El2  E22+  El4  C22+  El4  E22  (4.S9) 

A  similar  expression  must  be  obtained  for  other  combinations  of  values  of  k  and  1  as 
well  as  the  condition  that  the  RB  is  1.  This  produces  a  total  of  64  possibilities  of 
received  symbol  states  that  result  in  a  correct  RB.  A  computer  generated  table  was 
used  ,in  order  to  identify  the  equivalent  error  patterns  for  other  cases  involving  MB=  1 
and  LB=  1.  The  appropriate  listing  are  shown  in  Appendix  C.  Similai  analyses  must 
be  carried  out  in  order  to  determine  the  unconditional  probabilities  Pr(MB  correct)  as 
well  as  Pr(LB  correct}.  Such  analyses  demonstrate  that  due  to  the  inherent  symmetries 

Pr{  LB  correct  }  =  Pr{  MB  correct  }  =  Pr{  RB  correct  }  a  PE  (4.90) 

so  that  using  Equation  4.86  along  with  similar  forms  obtained  for  k  =  1  and  1=  1  yield 

PE  =  1  -  1:2  [  2C14  C22  +  9  Cl2  C22  El2+  Cl4  E22  +  3  Cl4  E2  C2^ 

+  10  Cl2  El2  E2:+  2  El4  C22  +  El4  E22  +  17  Cl2  C2  El2  E2  + 

+  Cl3  C22  El  +  3  Cl3  El  E22  +  2  Cl  El3  C22  +  2  Cl  El3  E22  + 

+  3  C2  El4  E2  +  4  Cl  C2  El3  E2  +  4  Cl3  C2  El  E2  ]  (4  91) 

In  Figures  4.4  and  4.5,  PE  (given  by  Equation  4.91  )  vs  SNR  have  been  plotted  for 
various  values  of  a  and  £  respectively. 

C.  NUMERICAL  RESULTS  AND  DISCUSSION  OF  RESULTS  FOR  8-PSK 
AND  8-DPSK  DBDR'S  WITH  RECEIVER  LOCAL  OSCILLATOR  PHASE 
ERROR. 

Previous  sections  of  this  chapter  have  focused  on  the  performance  of  8-PSK  and 
8-DPSK  DBDRs  operating  in  the  presence  of  a  receiver  local  oscillator  phase  error. 
Without  any  phase  error,  we  have  shown  the  performance  plots  of  S-PSK  and  S-DPSK 
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Figure  4.5  PE  Versus  SNR  of  3-DPSK  DBDR  for  Various  Values  of  e  (a 


DBDRs  on  Figure  4.3,  thus  demonstrating  the  performance  penalty  of  the  latter 
receiver.  While  usage  of  8-DPSK  provides  good  bandwidth  efficiency  (same  as  8-PSK) 
and  relatively  good  noise  immunity,  the  BER  of  8-PSK  is  lower  than  that  of  8-DPSK 
for  £~  0°,  as  expected 

The  value  of  using  8-DPSK  versus  8-PSK  can  be  appreciated  more  fully  when 
comparing  the  performance  plots  of  Figures  4.1  and  4.5.  These  show  that  if  the 
receiver  has  a  local  oscillator  phase  error  restricted  to  0^£^66.5°,  for  similar  SNR 
values,  the  BER  of  the  8-DPSK  receiver  is  lower  than  that  of  the  8-PSK  receiver. 
Since  such  phase  errors  can  be  typically  expected  to  be  small,  this  is  a  significant 
advantage.  However  it  must  be  pointed  out  that  for  phase  errors  of  more  than  10°, 
both  receivers  exhibit  severely  degraded  performance  that  would  make  both  systems 
unsuitable  for  operation  in  a  practical  system.  The  8-DPSK  receiver  does  not  exhibit 
complete  insensitivity  to  phase  errors  as  encountered  in  phase  measurement  receivers. 
Only  when  a  complete  phase  revesal  takes  place  (i.e.,  £=  180°)  is  the  performance  of 
the  8-DPSK  DBDR  unaffected,  while  the  8-PSK  receiver  produces  a  nearly  useless 
output. 

The  angle  a  set  to  22.5°  is  optimum  in  minimum  BER  sense  for  the  8-PSK 
DBDR.  The  variations  in  PE  versus  a  are  shown  in  Figure  4.6,  for  £=  0°  for  8-PSK, 
and  a  similar  plot  of  PE  versus  a  is  plotted  for  various  values  of  SNR  and  £=  0°  on 
Figure  4.7  for  the  8-DSPK  case. 


V.  PERFORMANCE  ANALYSES  OF  16-PSK  AND  16-DPSK  DBDR'S 
WITH  LOCAL  OSCILLATOR  PHASE  ERROR 

In  this  chapter,  DBDRs  will  be  analyzed  for  16-PSK  and  16-DPSK 
communication  systems  operating  in  the  presence  of  a  receiver  local  oscillator  phase 
error.  The  analysis  procedure  will  be  similar  to  that  carried  out  in  the  previous  chapter 
for  S-PSK  and  8-DPSK  systems. 

The  performance  of  16-PSK  DBDR's  has  been  investigated  and  the  methodology 
for  carrying  out  the  performance  analyses  have  been  well-documented  by  Bukofzer 
fRef.  8],  In  this  research,  these  analyses  will  be  extended  to  the  16-DPSK  system  as 
well  as  the  16-PSK  system  modified  to  operate  in  the  presence  of  a  receiver  local 
oscillator  phase  error. 

A.  16-PSK  DBDR  PERFORMANCE  ANALYSIS  IN  THE  PRESENCE  OF 

RECEIVER  LOCAL  OSCILLATOR  PHASE  ERROR 

The  possible  transmitted  signals,  which  are  represented  as  vectors  in  Figure  1.2.b. 

are  mathematically  given  by  Equation  4.1  with  i=  1,2,3 . 16.  The  same  notation 

introduced  in  Chapter  IV  will  be  used  to  represent  the  signals  using  the  orthonormal 
set  given  by  Equation  4.4.  Therefore  the  transmitted  signals  are  expressed  as 

Sj(t)=  <pj(t)  VE^  cos0j(t)  +  <p2(t)  sin0j(t)  (5.2) 

where  i=  1,2,.. .,16,  0:St^Ts  and 

/a  +  (i-1)  Jt/8  i=  1,5,9,13 
P  +  (i-2)  71/ 8  i=  2,6,10,14 

0.(t)  =  (i+1)  7t/s  -P  i=  3,7,1 1,15  (5.3) 

i  7t/8  -a  i  =  4,8,12,16 

\ 

with  a  and  P  as  shown  on  Figure  1.2.b 

Assuming  all  signals  are  equally  likely  to  be  transmitted,  for  the  16-PSK  case 
Equation  4.7  can  be  written  as 
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(5.4) 


Pr{  bit  correct 


16 

=  1/16  y  Pr(  bit  correct  /  Sj(  t)  transmitted  } 
i=  1 


The  notation  to  denote  the  four  bits  that  make  up  a  data  state  has  been  introduced  in 
Chapter  I.  Also  the  notation  introduced  by  Equation  4.13  is  used  here,  namely 

Pr{  bit  correct  /  s;(t) }  =  Pr{  bit  correct  /  s;(t)  transmitted  }  i=  1,2. .,16  (5.5) 


From  Figure  1.2.b,  the  conditional  probabilities  for  the  correct  reception  of  the  LB, 
MBL,  MBR,  RB  are  given  by 


Pr{  LB  correct  /  s;(t) }  = 


Pr{  r,  .  >  0  /  s.(t) }  i=  1,2,3,4,13,14,15,16 
1  ,1  1 


(5.6) 


( 

Pr{ 

rl,i<0  / 

s;(t) } 

i=  5,6,7,8,9,10,11,12 

>r{ 

r  2-r  2 
rl,i  r2,i 

V 

o 

—V5 

}  i- 1,2,7,8,9,10,15,16 

Pr{ 

MBL  correct 

/  Sj(t)  }  =< 

(5.7) 

Pr{ 

» 

r  2-r  2 
rl,i  r2,i 

A 

O 

—V* 

r* 

}  i=  3,4,5,6,11,12,13,14 

Pr{ 

VTh>0 

/  S;(t)  } 

i-  1,4,5,8,9,12,13,16 

Pr{ 

MBR  correct 

/  Sj(t)  }  = 

(5.8) 

Pr{ 

VTh<0 

/  Sj(t)  ) 

i=  2,3,6,7,10,1 1,14,15 

[Pr{ 

r2,i>0/ 

S;(t)  } 

i=  1,2, 3, 4, 5, 6, 7, 8 

Pr{ 

RB  correct  / 

S;(t)  }  =  ' 

(5.9) 

iPr{ 

r2,i  <  0  / 

s-/*) ) 

i  =  9, 10, 11, 12, 13, 14, 15, 16 

where  the  rk  ;  are  given  by  Equation  4.9  and 


Th 


2  \2 


=  <  ri,i  - V  > 


<rl,i  f2,i)2 


i- 1,2,3 . 16 


(5.10) 


The  analyses  methods  which  have  been  given  in  Chapter  IV,  can  be  applied  to 
determining  the  conditional  error  probabilities  for  LB,  MBL,  RB.  For  MBR,  the 
analysis  procedure  which  has  been  used  by  Bukofzer  [Ref.  8],  will  be  utilized.  In  order 
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l. 

•t 

t> 


to  define  the  conditional  probability  sets  of  Equation  5.8,  additional  steps  must  be 
completed  as  demonstrated  in  the  sequel.  If  s;(t)  transmitted,  the  YIBR  is  correctly 
recovered  if  for  i=  1,4,5,8,9,12,13,16,  Vj^>  0  and  if  for  i=  2,3,6,7,10,1 1,13,14,  Vy^  <  0. 
Even  thoueh  the  r  •  are  conditional  Gaussian  random  variables,  because  of  the  highly 
nonlinear  form  of  Equation  5.10,  a  more  productive  approach  to  the  determination  of 
the  conditional  p.d.f.  of  Vy^  has  been  introduced  (see  Reference  8  for  details)  rather 
than  using  transformations  of  Gaussian  random  variables.  By  introducing  the 
following  notation,  namely 


Tj  j=  component  ofr(t)  along  <p1(t)  =  £  cos  n 
r-,  .=  component  of  r(t)  along  <p-,(t)  =  £  sin  n 
it  can  be  demonstrated  that 

Vy^  =  £4  cos  4n 

where 

£  =  length  of  the  vector  =  J(r.  .)2  +  (r,  .)2 
q  =  vector  angle  that  with  respect  to  the  positive  x  axis 
=  arctan  (r2  •  I  r^) 


(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 


Since  £4>0,  it  is  clear  that  the  test  involving  Vy^  being  greater  than  or  less  than  zero 
is  now  equivalent  to  the  test  involving  cosdij  in  the  same  manner.  This  implies  that 
the  p.d.f.  of  n  must  be  obtained  which,  as  shown  in  Reference  8,  leads  to  integrals  that 
cannot  be  evaluated  in  closed  form  and  therefore  must  be  evaluated  via  numerical 
integration. 

By  utilizing  the  analysis  procedure  from  Reference  8  and  a  methodology 
introduced  in  the  previous  chapter  for  the  8-PSK  case,  the  conditional  probabilities  for 
the  correct  reception  of  a  bit  in  the  presence  of  receiver  local  oscillator  phase  error  is 
given  in  shorthand  form  by  the  following  expressions 


Pr{  LB  correct  /  s£( t)  }  =  1  -  Q{  Es  XQ  f  (a.p.c)  },  i=  1,2,3...  16 


(5.16) 


where 


f(a.p,c)  = 


cos  (a-e) 

i=l,9 

COS (P-E) 

1=2,10 

sin  (P  +  c) 

i  =  3,ll 

sin  (a  +  £) 

i  =  4,12 

sin  (P-e) 

i  =  5,13 

sin  (a-£) 

i  =  6,14 

cos (P  +  £) 

i  =  7,15 

cos  (a  +  £) 

i  =  8,16 

(5.17) 


and  a  and  p  are  shown  on  Figure  1.2.b.  Furthermore 

Pr(  MBL  correct  /  s;(t)  }  =  Q{-v/  Es/N0  (cos  y  +  sin  y)  .  Qf-v'Ej/Ng  (cos  y-sin  y)  }  + 

+  Q{VEs;X'0  (cos  y-sin  y) ).  Q(VES  XQ  (cos  y  + sin  y)  )  (5.18) 


where 


y-f(a-e)  i=  1,5,9,13 
l(P-€)  i=  2,6,10,14 


(5.19) 


Pr(  MBL  correct  /  Sj(t) }  = 


=  1  -  Q{-VES/N0  (cos  y  +  sin  y) .  Q{-v/Es,'N0  (cos  y-sin  y)  }  + 

+  Q{VES/N0  (cos  y  +  sin  y)}  .  Q{VES ,N0  (cos  y-sin  y )}  (5.20) 


where 


y  =  f  (P-€)  i=  3,7,1 1,15 
l(a-E)  i  =  4,8,12,16 


(5.21) 


For  the  calculations  involving  the  conditional  probability  for  the  MBR,  the  joint 
p.d.f.  of  the  random  variables  q  and  £  has  been  determined  by  Bukofzer  [Ref.  8] 
namelv 


re.n  s, 


(L.H)  - 


| L|  f_  r  (L  cosH  ,  L  sinH)  + 
l.i’  2,i;s. 

+  |L|  1-  r  '  (-L  cosH  ,  -L  sinH) 

ri,i’r2,i/s. 


(5.22) 


m  "  «  *  I  *  «  *  »  *  JH  *  q  *«  •  m  *  *  '  |  *  »  *  b  ^  j|^  k**  ■  *  a  a^*  a**  .*»  .N  *»  \ 


WWW 


which  must  be  integrated  to  produce  the  marginal  p.d.f.  of  rj-  This  results  in 

fn;s_  (H  '  s.)  =1.2  exp(-Es  N0)  +  {  1  +  V  4E$ ..JtN0  . 

*.  (  exp[  ES;N0  cos2  (H-£j)l  cos(H-^;) .  Q[ -Jl  Es  NQ  (cos(H-£.)  ]  } 


(5.23) 


where 


=  0j(t)  -  e  i=  1,2.... 16 


(5.24) 


where  9.(t)  is  given  by  Equation  5.3  Therefore,  the  probability  that  Vj^>0  involves 
computing  the  probability  that  cos4r|>0.  Observe  that  cos4q>0  for  -7t;  8  <  x\  <  n  8, 


3tc,  S  <  r|  <  5tc  S.  In.  8  <  r\  <  9n.  8,  and  1  Itt  8  <  <  13ti/S.  Therefore  the  calculation  of 
the  conditional  probability  of  the  MBR  being  correct  given  that  Sj(t)  was  transmitted 
can  be  calculated  by  integrating  the  conditional  p.d.f  of  rj  over  the  appropriate 
regions,  namely 

Pr(  MBR  correct  /  s^t)  }  =  Pr{  cos4n>0  /  Sj(t)  }  = 
tt/8  5rc;8  9tt/8 

-  1  Vs, lH  si]  dH  +  1  fn;s1  lH  sii dH  + 1  Vs,  |H  si]  dH  + 

-n:  8  3)1/8  77C/8 

1371/8 

+  J  fn/sj  ^H;siJ dH 

llTt  8 

7t.  8 

=  j  {  fn/  [H  /  Sj]  +  fn  $1  [(H  +  Tt/2)  /  Sj]  +  fn/s  [(H  +  ti)  /  Sj]  +  (5.25) 

-71/ 8 

+  fn/  [(H  + 371/2)  /  Sj] )  dH 

By  utilizing  the  appropriate  form  for  the  conditional  p.d.f.  of  x\  when  Sj(t)  is  assumed 
transmitted,  results  in 

Pr(  MBR  correct  s.(t)  }  =  1  2  exp(-Es  NQ)  + 

71  8 

+  v'  Es  TtNg  j  (  exp{-  E$  Nq  sin-(H-X)]  cos(H-X)  .  {Q[-v2  Es  NQ  (cos(H-X)]  - 

-71  8 
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-  Q[v/2  Es, N0  (cos(H  +  X.)] }  +  (  exp[-Es, \Q  cos2(H-X)]  sin(H-X)  Qf-V2  Es  NQ  (sin(H-X)j  - 
-Q[V2  Es/N0(sin(H +  *.)]}  }  dH  (5.26) 

where 


a  +  £ 

i=  4,8,12,16 

a-£ 

i=  1,5,9,13 

(5.27) 

P-E 

i==  2,6,10, 14 

P  +  £ 

i=  3,7,11,15 

Finally,  it  can  be  seen  that 

Pr{  RB  correct  /  Sj(t)  }  =  Pr{  LB  correct s;(t)  )  (5.28) 

so  that  by  using  Equations  5.4  and  5.16-5.28  for  correct  reception  probabilities  of  the 
individual  bits,  the  overall  BER  can  be  obtained  by  using 

Pr{  bit  in  error  }  =  1  -  Pr{  bit  correct  }  (5.29) 

This  results  in  the  (unconditional)  probabilities 

Pr{  LB  in  error  }  =  Pr{  RB  in  error  } 

=  1,  8  {  Q[  -J2  Es/N0  cos(a-e)]  +  Q[  J2  Es:N0  cos( P-e)]  + 

+  Q[  V2  ES'N0  sin(P  +  £)]  +  Q(V2  Es/N0  sin(a  +  £)l  + 

+  Q[  V^/Ng  sin(p-E)]  +  QtvTETNg  cos(P  +  £)1  + 

+  Q  (v^  Ej  Ng  cos(a  +  £)]  +  Q  [V2  Es  Ng  sin(a  +  £)]}  (5.30) 


Pr{  MBL  in  error  }  =  1/2  +  1/4  [  Q  (-VEs/N0  [cos(a-£)  +  sin(a-E)])  + 

+  Q  {-  V’Ej/Ng  [cos(a  +  £)  -  sin(a  +  £)|}  +  Q  {-VE^/Ng  [cos(P-£)  +  sin{ P-£)]}  + 
+  Q  {-v'Ej./Ng  [cos(P-£)  -  sin(P-£)|)  -  Q  (-v'Es  NQ  (cos(P  +  £)  +  sin<P  +  £)|}  - 
-  Q  {-VES/N0  (cos(p  +  £)  -  sin(p  +  £)|}  -  Q  {-V'ES,  NQ  [cos(a  +  £)  +  sin(a  +  £»}  - 


-  Q  {-VEj.  Ng  (cos(a  +  £)  -  sin(a  +  £)]}  ]  +  1,2  [  Q  {-VE~N0  [cos(a  +  £)  +  sin(a  +  £)]} 
.  Q  {-v'Es.N0  [cos(a-t)  -  sin(a-E)]}  +  Q  {-v'Es,N0  [cos(p  +  E)  +  stn(P  +  £)]}  . 

.  Q  r-s/Ej  N0  [cos(P  +  £)  -  sin(P  +  £)]}  -  Q  {-v'Es,N'0  [cos(a-£)  +  sin(a-£)]}  . 

•  Q  (-vEs;N0  [cos(a-£)  -  sin(a-£)]}  -  Q  {VES.N„  [cos(P-e)  +  sin(p-E)]}  . 

•  Q  {-VEg/Ng  [cos(p-£)  -  sin(p-£)]}  ]  (5.31) 

For  the  MBR  the  result  appears  in  integral  form.  That  is 


Pr{  MBR  in  error }  =  1  -  1/2  J  Es/7tN'0  J  {  exp{-Es/N'0  sin2[H-(a-£)j}  cos[H-(a-£)]. 


-  Q[-VEs/N0  cos(H-(a-£))]}  +  exp{-E$/N0  cos2[H-(a-£)]}  sin[H-(a-£)] 

-  QI-VEj/Ng  sin(H-{a-£))]}  +  exp{-Es  \'Q  sin:[H-(p-£)])  cos[H-(p-£)] 

-  Q[-VES/N0  cos{H-(p-£))]}  +  exp{-Es;N0  cos2[H-(p-£)]}  sin[H-(P-e>] 

-  Q[-VEj/N0  sin(H-(p-£))]}  -  exp{-Es;N0  sin2[H  +  (p  +  £)]}  cos[H  +  (p  +  £)] 

-  QI-VEj/Nq  cos(H  +  (P  +  £))]}  -  exp(-Es'N0  cos2[H  +  (P  +  e)]}  sin[H  +  (P  +  £)] 

-  Q[*VEs/N0  sin(H  + (a +  £))]}  +  exp(-Es  NQ  sin2[H  +  (a  +  £)]}  cos[H  +  (a  +  e)] 

-  Q[-V/Es,  Nq  cos(H +  («  +  £)]}  +  exp(-Es.  N0  cos2(H  +  (a  +  £)]}  sin[H  +  (a  +  £)] 

-  QI-TeT^  sin(H  +  (a  + 1»]>  }  dH  (5.32) 


The  overall  BER  can  now  be  determined  from  Equations  5.30,5.31  and  5.32  by 
evaluating  the  final  result  for  PE,  the  bit  error  probability,  namely 


PE  =  1/4  (  Pr{  LB  in  error  }  +  Pr{  MBL  in  error  } 
+  Pr{  MBR  in  error  }  +  Pr{  RB  in  error  }  } 


(5.33) 


It  can  be  seen  that  Equation  5.33  consists  of  many  terms  that  by  themselves  do  not 
yield  as  much  insight  as  the  computed  values  of  PE  presented  in  plotted  form. 
Therefore  plots  of  BER  vs.  SNR  for  various  values  of  a  and  £,  are  given  as  Figures  5.1 
and  5.2.  While  a  and  P  can  be  set  independently,  from  an  optimization  point  of  view, 
a  +  p  =  7r  4  has  been  shown  to  produce  the  lowest  BER  for  given  values  of  SNR  when 
£  =  0°  [Ref.  S] 
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Figure  5.1  PE  Versus  SNR  for  16-PSK  DBDR  (o-  11.25°,  c*0°). 
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B.  16-DPSK  DBDR  PERFORMANCE  ANALYSIS  IN  THE  PRESENCE  OF 
RECEIVER  LOCAL  OSCILLATOR  PHASE  ERROR 

Similar  analyses  on  DBDR  performance  have  been  carried  out  for  16-DPSK  (see 
Chapter  IV. B),  so  that  the  notation  introduced  previously  will  be  used  here.  That  is. 

Ck  =  Pr  (  k  correct  }  ,  Ek  =  Pr  (  k  in  error  }  (5.34) 

where  now  k=  1  for  LB,  k=  2  for  MBL,  k=  3  for  MBR  and  k  =  4  for  RB. 

In  order  to  evaluate  Pr  {  LB  correct  }  with  differential  decoding,  we  have  to 
evaluate 

Pr{  LB  correct )  =  1:2  [  Pr{  LB  correct  /  LB  =  0  }  + 

+  {  Pr{  LB  correct .  LB=  1  }  j  (5.35) 

where 

Pr{  LB  correct  /  LB  =  0  }  =  (5.36) 

1  1  1 

=  n  I  Pr{  LB  correct/ LB  =  0,  MBL  =  i,  MBR  =  j,  RB  =  k}  x 
i  =  0  j  =  0  k  =  0 

x  Pr{  MBL  =  i,  MBR  =  j,  RB  =  k} 

where  for  i,  j,  k=  0,  1 

Pr{  MBL=i,  MBR  =  j,  RB  =  k  )  =  1/8  (5.37) 

In  Chapter  IV,  it  has  been  explained  that  with  DPSK  modulation  in  conjuction 
with  bit  direct  detection  methods,  certain  bit  errors  in  the  past  and  presently  recovered 
symbol  states  do  not  result  in  errors  in  the  recovered  data  states.  Such  bit  errors  in 
successive  symbol  states  must  be  considered,  because,  even  though  symbol  errors  may 
occur,  the  delivered  data  may  still  be  correct.  Therefore,  these  errors  have  to  be 
accounted  for  and  considered  in  the  computation  of  the  probability  of  a  bit  being 
correct.  Figure  5.3  shows  one  example  of  the  situation  being  discussed.  That  is,  the 
data  state  (0,0, 0,0)  is  assumed  to  occur  and  the  symbol  state  (0,0, 0,0)  is  assumed  to 
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Figure  5.3  Example  of  Error  Patterns  for  16-DPSK  System  Where  Data  State. 
Previous  Transmission,  and  Present  Transmission  are  0000  with  RB  Correctly  Recovered. 
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Figure  5.3  Continued. 
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have  been  previously  transmitted.  This  means  that  according  to  the  encoder  matrix 
shown  in  Chapter  II,  the  symbol  state  (0.0, 0,0)  must  now  be  generated  for 
transmission.  Figure  5.3  shows  the  receptions  of  the  symbol  states  which  produce 
recovered  data  states  in  which  the  RB  is  correctly  delivered.  Observe  that  there  are  12S 
such  combinations  which  yield  0  for  the  RB.  Therefore  notation  has  to  be  introduced 
to  account  for  error  patterns  of  the  above  type,  namely  Pr{  RB  correct,  (g,h,p,q)  / 

LB  =  m,MBL  =  i,MBR  =  j,RB  =  k.  },  where  g,h,p,q=  0,1.  The  meaning  of  this 
expression  is  the  probability  of  the  RB  being  correctly  recovered  when  the  data  state  is 
(m,j,i,k)  and  the  delivered  data  state  is  (g,h.p,q)  where  q  =  0  in  this  case  (otherwise  the 
RB  could  not  be  correct  due  to  the  transmission  of  the  RB  =  0).  Therefore 

Pr{  RB  correct  /  LB  =  0,  MBL  =  0,  MBR  =  0,  RB  =  0  }  = 

1  1  1 

Z  Z  Pr{  RB  correct,  (g,h,p,0) ;  LB  =  0,  MBL  =  0,  MBR  =  0,  RB  =  0  }  (5.3S) 

g  =  0  h  =  0  p  =  0 

Since  each  element  of  the  symbol  state  is  independently  recovered,  we  see  from  Figure 
5.3  that 

Pr{  RB  correct, (0,0, 0,0)  /  LB  =  0,MBL  =  0.MBR  =  0,RB  =  0  }  =  Cl2  C22  C32  C42  + 

+  Cl2  C22  E32  C42  +  Cl2  C22  E32  E42+  Cl2  E22  C32  E42+  Cl2  E22  E32  C42  + 

+  Cl2  E22  E32  E42+  El2  +  C22  C32  C42+  El2  C22  C32  E42-i-  El2  C22  E32  C42  + 

+  El2  C22  E32  E42+  El2  E22  C32  C42+  El2  E22  E32  C42+  El2  E22  E32  E42  + 

+  Cl2  C22  C32  E42  (5.39) 

However,  from  Equations  5.34,  and  5.30,  it  is  shown  that  Cl=  C4  and  El=  E4. 
Therefore, 

Pr{  RB  correct,  (0,0, 0,0)  /  LB  =  0,  MBL  =  0,  MBR  =  0,  RB  =  0  )  = 

=  Cl4  C22  C32+  2  Cl2  C22  C32  El2  +  C22  C32  El2  E22 
+  2  Cl4  C22  E32  +  2  Cl2  C22  El2  E32+  Cl4  C32  E22  +  El4  E22  E32 
+  2  Cl2  C32  El2  E22  +  C22  El2  E22  E32  +  Cl2  El2  E22  E32  + 

+  C32  El4  E22  +  Cl2  C32  El4  (5.40) 
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A  similar  expression  must  be  obtained  for  other  combinations  of  values  of  g,h  and  p  as 
well  as  the  condition  that  RB  is  1.  This  produces  a  total  of  256  possibilities  of  received 
symbol  states  that  result  in  a  correct  RB.  Computer  generated  outputs  were  used  to 
identify  the  equivalent  error  patterns  for  other  cases  involving  LB=  1,  MBL=  1  and 
MBR  =  1.  These  results  are  shown  in  Appendix  C.  Similar  analyses  must  be  carried 
out  in  order  to  determine  the  unconditional  probabilities  Pr(LB  correct},  Pr(MBL 
correct}  as  well  as  Pr{MBR  correct}.  Such  analyses  demonstrate  due  to  inherent 
symmetries  that 

Pr{LB  correct}  =  Pr{MBL  correct}  =  Pr(MBR  correct}  =  Pr{RB  correct}  =  PE 
so  that  ubing  Equation  5.38  along  with  similar  forms  obtained  for  k=  1  and  1=  1  yield 

PE  =  1  -  1/2  {  2  Cl4  C22  C32  4-  Cl4  C22  C3  E3  +  3  Cl4  C22  C3  E2  + 

+  2  Cl4  C2  C3  E2  E3  +  6  Cl3  C22  C32  El  +  4  Cl3  C22  C3  El  E3  + 

+  12  Cl3  C2  C32  C3  El  E2  +  8  Cl3  C2  C3  C3  El  E2  +  6  Cl3  C22  El  E32  + 

+  6  Cl3  C22  C32  El  +  Cl4  C3  E22  E3  +  4  Cl3  C3  El  E22  E3  + 

+  2  Cl4  El2  E32  +  Cl4  C22  E32  4-  3  Cl4  C2  E2  E32  +  2  C22  C32  El4  + 

+  Cl4  C32  E22  +  10  Cl2  C22  C32  El2  +  10  Cl2  C22  El2  E32  + 

+  10  Cl2  C32  El2  E22  4-  8  Cl2  El2  E22  E32  4-  18  Cl2  C2  C32  El2  E224- 
+  18  Cl2  C2  El2  E2  E32  4-  6  Cl3  El  E22  E32  4-  12  Cl3  C2  El  E2  E32  + 

4-  12  Cl2  C2  C3  El2  E2  E3  4-  6  Cl2  C22  C3  El2  E3  +  5  Cl2  C3  El2  E22  E3  4- 
4-  6  Cl  C22  C32  El3  +  6  Cl  C32  El3  E22  +  6  Cl  C22  El3  E22+  C22  El4  E32 
4-  12  Cl  C2  C32  El3  E2  4-  12  Cl  C2  El3  E2  E32  +  6  Cl  El3  E2,  E32  + 

+  8  Cl  C2  C3  El3  E2  E3  4-  4  Cl  C3  El3  E22  E3  +  4  Cl  C22  C3  El3  E3  + 

4-  C32  El4  E224-  2  El4  E22  E32  4-  C22  C3  El4  E3  4-  C3  El4  E22  E34- 
4-  3  C2  El4  E2  E32  +  2  C2  C3  El4  E2  E3  +  2  C2  C32  El4  E2  }  (5.42) 

In  Figure  5.4,  PE  (given  by  Equation  5.42)  is  plotted  for  various  values  of  a.  Also 

Figure  5.5  includes  BER  for  various  values  of  e  when  a=  11.25° 
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C.  NUMERICAL  RESULTS  AND  DISCUSSION  OF  RESULTS  FOR  16-PSK 
AND  16-DPSK  DBDR'S  WITH  RECEIVER  LOCAL  OSCILLATOR  PHASE 
ERROR 

In  this  chapter  16-PSK  and  16-DPSK  DBDRs  have  been  analyzed  when 
operating  in  the  presence  of  a  receiver  local  oscillator  phase  error.  The  performance 
results  which  have  been  obtained  for  the  16-ary  case,  are  similar  to  those  obtain  for 
the  8-ary  case.  Hence,  instead  of  repeating  the  conclusions  similar  to  these  presented 
for  the  8-ary  case,  conclusions  that  are  unique  to  the  16-ary  case  will  be  presented  in 
this  part.  Main  differences  between  8-ary  and  16-ary  cases  which  are  observed  via 
analyses  and  numerical  evaluations  are 

1.  The  optimum  angles  in  minimum  BER  sense,  a  and  P  between  the  symbols 
becomes  11.25°  and  33.75°  respectively  for  16-PSK  case,  as  well  as  for 
16-DSPK  case,  with  a  +  P  =  45° 

2.  The  usage  of  16-PSK  provides  better  performance  than  16-DPSK  for  £<33.5°, 
however  the  actual  BER  values  are  too  high  for  use  in  a  practical  system  when 
£  greater  than  a  few  degrees. 

3.  The  16-DPSK  system  performs  better  than  the  16-PSK  system  for  £>33.5°. 
However  again  the  actual  values  of  PE  under  these  circumstances  is  too  high 
for  a  practical  system. 

4.  The  16-DPSK  system  is  insensitive  to  phase  errors  of  7t  radians. 

Performance  comparisons  between  the  16-PSK  and  16-DPSK  DBDRs  are  plotted 

on  Figure  5.6,  and  BERs  are  plotted  on  Figures  5.7  and  5.8  respectively  for  various 
values  of  a  with  £=  0° 
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VI.  CONCLUSIONS 

In  this  research  the  use  of  DBDRs  had  been  analyzed  for  M-PSK  and  MDPSK 
communication  systems.  The  main  assumptions  in  the  analyses  are  a  Gray  coding 
mapping  for  the  symbols  and  an  Additive  White  Gaussian  Noise  channel  model.  It  is 
further  assumed  that  the  receiver  has  a  known  phase  error  in  its  local  oscillator  and  its 
effect  on  receiver  performance  is  analyzed  for  the  the  special  cases  in  which  M  =  8  and 
\1  *  16. 

In  addition  to  the  analytical  results,  graphical  results  are  presented  in  the  form  of 
plots  of  va.  luLons  in  PE  for  various  values  of  phase  angle  a,  SNR.  and  local  oscillator 
phase  error  z. 

Optimum  conventional  M-PSK  and  M-DPSK  receivers  are  designed  to  detect 
symbols  while  DBDRs  are  designed  to  detect  bits.  However,  local  oscillator  phase 
errors  degrade  the  performance  of  M-PSK  and  M-DPSK  conventional  receivers.  Part 
of  the  work  undertaken  was  intended  to  investigate  the  effect  of  such  phase  errors  on 
DBDR's.  The  variations  of  PE  versus  £  are  plotted  in  Figure  6.1  for  8-PSK  and 
S-DPSK,  and  in  Figure  6.2  for  16-PSK  and  16-DPSK.  If  £  (receiver  local  oscillator 
phase  error)  can  be  constrained  to  a  small  interval,  namely  0°<£<30°,  the  M-PSK 
system  provides  better  performance  than  the  M-DPSK  system.  For  values  of  E  outside 
this  range,  both  systems  have  high  error  rates  that  make  the  receivers  unsuitable  for 
practical  applications.  However,  for  complete  phase  reversal,  i.e.,  £=  180°,  the  M- 
DPSK  system  suffers  no  performance  penalty  whereas  the  M-PSK  system  can  not 
properly  operate. 

The  performance  degradation  due  to  a  phase  error  has  been  analyzed  also  from 
a  different  point  of  view.  That  is,  when  E  *  0°,  there  is  an  increase  in  PE  for  both  the 
M-PSK  and  M-DPSK  receivers.  The  ratio  of  the  phase  error  is  given  by 


PER  =  PE0  PEj 


where 


PEq=  Probability  of  bit  error  @  £  =  0° 
PEj  =  Probability  of  bit  error  g  £  *  0° 


& 
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PER  has  been  analyzed  under  various  conditions  for  M-PSK  and  M-DPSK.  The 
results  have  been  plotted  in  Figures  6.3and  6.5  for  M  =  8,  and  in  Figures  6.6-6.S  for 
M=16,  as  a  function  of  SNR  for  different  values  of  £.  The  plots  demonstrate 
significant  differences  in  the  performance  degradation  of  the  M  =  8  versus  M  =  16  case. 
However,  for  high  SNR  values,  these  differences  tend  to  vanish. 

The  analysis  results  lead  to  the  following  conclusions 

1.  By  using  DBDRs,  the  error  probability  of  each  bit  in  the  symbol  state  can  be 
calculated  independently  and  more  importantly  these  error  probabilities  can  be 
varied  by  changing  the  angle  between  symbols. 

2.  For  optimum  BER,  a  should  be  set  to  22.5°  for  8-PSK  and  8-DPSK  systems, 
while  16-PSK  and  16-DPSK  systems  require  a=  11.25°  and  p  +  a  =  45°. 

3.  The  solution  for  the  phase  ambiguity  problem  cannot  be  achieved  by  using  M- 
DPSK  DBDR's,  unless  the  phase  ambiguity  is  180°.  Therefore,  if  the  perfect 
coherency  conditions  cannot  be  achieved,  the  M-PSK  system  should  be  used 
whenever  DBDR's  are  used. 

4.  When  £  =  0°,  the  DPSK  system  exhibits  good  performance.  However,  if  the 
£—  0°  condition  can  be  guaranteed  there  is  no  need  to  use  the  differential 
encoding  scheme. 
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LEGEND 
□  8-DPSK 
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_ SNR  (DB) _ 

Figure  6.3  PER  (c  =  0o  e=5°)  Versus  SNR  for  S-PSK 
and  S-DPSK  DBDR's  (a=  22.5°). 
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_ SNR  (DB) _ 

Figure  6.4  PER.  (£  =  0°  c=  22.5°)  Versus  SNR  for  S-PSK 
and  8-DPSK  DBDR's  («=  22.5°). 
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Figure  6.5  PER  (£  =  0°:c  =  45°)  Versus  SNR  for  S-PSK 
and  8-DPSK  DBDR's  («  =  22.5°). 
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Figure  6.o  PER  ic  =  0°  e  =  5°)  Versus  SNR  for  16-PSK 
«nd  lo-DPSK  DBDR's  (a-  11.25°). 


Figure  6.7  PER  (£  =  0°  e  =  11.25°)  Versus  SNR  for  16-PSK 
and  16-DPSK  DBDR's  (a  =  11.25°). 


Figure  6.8  PER  (c  =  0°  c  =  45°)  Versus  SNR  for  16-PSK 
and  16-DPSK  DBDR's  («=  11.25°). 
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APPENDIX  A 

DBDR  INPUT/OUTPUT  GENERATION  PROGRAM  (BELIZ) 

I.  USER  GUIDE 

This  program  generates  all  the  possible  input  /  output  combinations  in  an  M- 
DPSK  communication  system.  The  program  runs  interactively  and  sends  outputs  to 
the  files 

1.  FILE  FT08F001 

2.  FILE  RCDATA 

3.  FILE  XMDATA 

4.  FILE  R 

5.  FILE  ERR. 

These  files  are  used  for  information  transfer  between  subroutines  but  are  not  required 
for  program  execution.  It  is  possible  to  generate  input  /  output  combinations  for 
8-DPSK  and  16-DPSK  by  using  this  program,  while  for  32-DPSK  or  64-DPSK 
systems,  dimensions  and  format  statements  must  be  changed  accordingly,  and  a  larger 
character  variable  set  in  the  subroutine  XMTR  must  be  created.  The  program  accepts 
only  the  integer  equivalents  of  your  data  or  symbol  states.  An  example  run  of  the 
program  is  shown  below 

run  beliz 

EXECUTION  BEGINS... 

PLEASE  SELECT  ONE  OF  THEM 

(1)  IF  YOU  HAVE  TRANSMITTER 

(2)  IF  YOU  HAVE  RECEIVER 

1 

PLEASE  SELECT  ONE  OF  THEM 

(1)  IF  YOUR  SYSTEM  8-DPSK 

(2)  IF  YOUR  SYSTEM  16-DPSK 

1 

ENTER  THE  INTEGER  EQUIVALENTS  OF  SYMBOL  STATES 
IN  THE  ORDER  OF  COUNTER  CLOCK  WISE 
ENTER  THE  STATE  #  1 

7 

ENTER  THE  STATE  #  2 

5 

ENTER  THE  STATE  #  3 

1 

ENTER  THE  STATE  #  4 
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3 


ENTER  THE  STATE  #  5 


2 

ENTER  THE  STATE  #  6 


0 

ENTER  THE  STATE  #  7 


4 

ENTER  THE  STATE  #  8 


6 

ENTER  THE  INTEGER  EQUIVALENTS  OF  THE 
PHASE  ANGLES  OF  YOUR  CODING  SYSTEM 
1*ALPHA=  ? 


0 

1 

3 

2 

6 

7 

5 


2*ALPHA=  ? 

3*ALPHA=  ? 
4 * ALPHA =  ? 

5*ALPHA=  ? 
6*APLHA=  ? 
7*ALPHA=  ? 
8*ALPHA=  ? 


4 

PLEASE  ENTER  THE  SELECTION  CORRECTLY 

ENTER  1  IF  YOU  NEED  BINARY  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  2  IF  YOU  NEED  DECIMAL  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  3  IF  YOU  WANT  TO  TURN  BACK  MAIN  PROGRAM. 


1 

PLEASE  ENTER  THE  SELECTION  CORRECTLY 

ENTER  1  IF  YOU  NEED  BINARY  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  2  IF  YOU  NEED  DECIMAL  FORM  OF  SYMBOL  AND 
DATA  STATES  INCLUDING  SYMBOL  OUTPUTS 

ENTER  3  IF  YOU  WANT  TO  TURN  BACK  MAIN  PROGRAM. 


3 


ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  FT08F001 


non  no  n n  nnnnnnnnnnnnnnnnnnnnnnn 


2.  SOURCE  CODE 


THIS  PROGRAM  SIMULATES  A  M-DPSK  COMMUNICATION  SYSTEM 
PROGRAM  RUNS  TOTATLY  INTERACTIVE  THEN  SENDS  TO  OUTPUTS  FILES; 

FILE  FT08F001 , FILE  RCDATA, FILE  XMDATA, FILE  R, FILE  ERR. 

THESE  FILES  ARE  USED  FOR  INFORMATION  TRANSFER  BETWEEN 
SUBROUTINES  BUT  ARE  NOT  REQUIRED  FOR  PROGRAM  EXECUTION. 

YOU  CAN  SIMULATE  8-DPSK  AND  16-DPSK  BY  USING  THIS  PROGRAM 
FOR  32-DPSK  OR  64-DPSK  SYSTEMS  YOU  HAVE  TO  CHANGE  DIMENSIONS 
AND  FORMAT  STATEMENTS  ACCORDINGLY, ALSO  YOU  HAVE  TO  EXPEND 
CHARACTER  VARIABLE  SET  IN  THE  SUBROUTINE  XMTR. 

OPERATION  SPECIFICATIONS . ; 

THE  SYMBOL  STATES  REQUESTEd  DEPEND  ON  YOUR  SYSTEM  AND  CAN 
BE  ANY  SEQUENCE.  MAIN  THING  IS;  YOU  HAVE  TO  SELECT  YOUR  REFERENCE 
THEN  ENTER  THEM  BY  ONE  STEP  INCREMENTS , ALSO  SAME  SPECIFICATION 
EXISTS  FOR  YOUR  ANGLES . PROGRAM  ONLY  EXCEPTS  THE  INTEGER  EQUIVALENTS 
OF  YOUR  STATE  OR  ANGLE  CODES. 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


VARIABLE  DECLERATION 

INTEGER  N, L, ORDER , TM(64 , 64) ,PS(64) , TRDATA(256) ,M,KK,K, THETA (64) 
INTEGER  MAT ( 64 ) ,KL,RMAT(64,64) ,TR(64) ,TH(64,64) , RCVMAT (64 , 64) ,NN 
INTEGER  QQU , QU , SQU , LD , I S , NT ( 1 , 8 ) . KM ( 1 , 8 ) , KN ( 1 , 8 ) 

INTEGER  RCV(64,64) , T(256) ,TD(256) ,NEXTST(256) ,PRESTA(256) ,FL 
INTEGER  A.B,Y,KD,KY(1,8) ,KH(1,8) ,KJ(1,8) , KB (1,8) 

CHARACTER^  DATA (64, 64) 

CHARACTER* 16  SET 


***  CALLING  SUBROUTINES  *** 

CALL  INFO  (ORDER, PS, THETA, NN,N,FL) 

IF  ( FL . EQ . 1 ) THEN 
CALL  XMTR  ( ORDER, FL) 

CALL  XMSTA  ( THETA, TM, PS , ORDER) 

CALL  XMDAT (TM , ORDER , THETA , NEXTST , PS , T , TD ) 

CALL  CONVER  (N, NEXTST, NN, ORDER, FL) 

WRITE (6 , 1 ) 

1  FORMAT (T9 , 1  ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  FT08F001 ' ) 
ELSE 

CALL  RCVR( ORDER) 

CALL  RCVSTA  (PS , ORDER , THETA, RCVMAT) 

CALL  RCVDAT ( TR , PS , ORDER , NEXTST , RCVMAT ) 

CALL  CONVER  (N, NEXTST, NN, ORDER, FL) 

CALL  ERROR ( FL , THETA , RCVMAT , PS , ORDER , TR , TM , T , NEXTST , TD ) 

WRITE (6, 2) 

2  FORMAT (T9, '  ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  FT08F001',/, 
*T7 , 1  ERROR  ANALYSES  RESULTS  WAS  SENT  TO  THE  FILE  ERR 1 ) 

END  IF 

STOP 

END 


SUBROUTINE  INFO ( ORDER , PS , THETA , NN , N , FL ) 

THIS  SUBROUTINE  PROVIDES  INFORMATIONS  ABOUT  OUR  SYSTEM 

INTEGER  N. L, ORDER, TM(64, 64) ,PS(64) , THETA (64) ,FL 
CHARACTER* 1  DATA(64,64) 

WRITE  (6,*)  1  PLEASE  SELECT  ONE  OF  THEM' 

WRITE  (6,*)  '(1)  IF  YOU  HAVE  TRANSMITTER' 

WRITE  (6  *)  '(2)  IF  YOU  HAVE  RECEIVER' 

READ  (6,*)  FL 

WRITE  (6,*)  1  PLEASE  SELECT  ONE  OF  THEM' 

WRITE  (6,*)  '(1)  IF  YOUR  SYSTEM  8-DPSK1 


nnnnn  o  nnonnn 


WRITE  (6  *)  '(2)  IF  YOUR  SYSTEM  16-DPSK' 

READ  (6,*)  L 

IF(L.EO.l)  THEN 

0RDER=8 

NN=63 

N=6 

END  IF 

IF(L.EQ.2)  THEN 

0RDER=16 

N=8 

NN=255 

ENDIF 

WRITE  (6,*)  ‘ENTER  THE  INTEGER  EQUIVALENTS  OF  SYMBOL  STATES' 
WRITE  (6,*)  'IN  THE  ORDER  OF  COUNTER  CLOCK  WISE1 
DO  11  K  =1. ORDER 
WRITE (6 ,20)  K 

20  FORMAT(l OK, 'ENTER  THE  STATE  #',I3) 

READ ( 6  ,  * )  PS (K) 

11  CONTINUE 

WRITE  (6,*)  'ENTER  THE  INTEGER  EQUIVALENTS  OF  THE  ' 

WRITE  (6,*)  'PHASE  ANGLES  OF  YOUR  CODING  SYSTEM  ' 

DO  33  L=l, ORDER 
WRITE  (6,10)  L 

10  FORMAT  ( 1QX , 13 , ' *ALPHA=  ?') 

READ ( 6 , * )  THETA ( L ) 

33  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  XMTR (ORDER , FL) 


THIS  SUBROUTINE  PROVIDES  GENERAL  TRANSMITER  MATRIX  FOR 
M-DPSK  SYSTEM. 

ITS  OUTPUT  CAN  BE  SEEN  ON  EITHER  FILE  FT08F001  OR  FILE 
XMDATA . 

INTEGER  N (L, ORDER, FL 

CHARACTER* 16  SET 

SET (1:16)=  ' ABCDEFGHI JKLMNOP ' 

OPEN (UNIT=11,FILE=' XMDATA' , STATUS= ' OLD ' ) 

IF  (FL.EQ.l)  THEN 
WRITE  (8, 81) ORDER 

81  FORMAT(///,T15, 'TRANSMITTER'S  ENCODER  MATRIX  ANALYSIS 

*  FOR1 ,13, ' -DPSK' 

* , / , 60 ( ' - ' ),//,T5, ' GENERAL  TRANSMITTER  MATRIX' ,/ ,T2,31( '-’) ) 

FT  QF 

WRITE  (8,88) 

88  FORMATS///, T15, 'TRANSMITTER  INFORMATION  FOR  THE  GIVEN  RECEIVER' 
*,/,72('-' ) ,//,T5, 'GENERAL  TRANSMITTER  MATRIX' ,/,T2,31( '-') ) 
ENDIF 

DO  3  N=l, ORDER 
IF  (N.EQ.l)  THEN 
WRITE  (§,*)  SET(N:ORDER) 

WRITE  (11,*)  SET(N:ORDER) 

ELSE 

WRITE  (8,*)  SET (N: ORDER) ,SET(1 :N-1) 

WRITE  (11,*)  SET(N:ORDER  ,SET(1:N-1) 

ENDIF 

3  CONTINUE 

CLOSE  (11) 

RETURN 

END 

SUBROUTINE  XMSTA ( THETA , TM , PS , ORDER ) 


THIS  SUBROUTINE  CONVERTS  GENERAL  XMTR  MATRIX  IN  TO  OUR  SPECIAL 
XMTR  SYMBOL  MATRIX. 

GENERAL  XMTR  MATRIX  IS  READ  FROM  FILE  XMDATA. 


INTEGER  TM(64,64) ,PS(64) , THETA ( 64 ), ORDER 

90 
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CHARACTER*!  DATA (64, 64) 

0PEN(UNIT=11 , FILE= 1 XMDATA1 , STATUS  ='OLD' ) 

IF  ( ORDER, EQ. 8) THEN 

WRITE (3 , 99 ) (THETA ( I ) ,1=1, ORDER) 

FORMAT (// / ,T5 , '  TRANSMITTER  MATRIX  '  ,  /  ,  31 ( 1  -  1  ) , // , 

*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS)' 
*,//,4X,3<I3),/,35('-')) 

READ  (11,10)  ( (DATA(L ,N) ,N=1 , 8) ,L=1 ,3) 

FORMAT  (IX, 3A1 ) 

ELSE 

WR I TE  ( 8 , 43 )  ( THETA (I), 1=1, ORDER) 

F0RMAT(///,T5, 'TRANSMITTER  MATRIX  1  ,  /  ,  31 (  1  -  1 ) , // , 

*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS)' 

*, //.4X,  16(13)  /, 700-')) 

READ  (11,20)  ( (DATA(L,N) ,N=1 , 16) ,L=1 ,16) 

FORMAT  (IX, 16A1) 

END  IF 

DO  3  K=l, ORDER 
DO  4  N=l, ORDER 
DO  6  M=  0, ORDER 

IF ( I CHAR ( DATA ( K , N ) ) . EQ . I CHAR ( DATA ( 1 , ORDER-M ) ) ) THEN 
TM(K,N)=PS(ORDER-M) 

ENDIF 

CONTINUE 

CONTINUE 

IF  ( ORDER. EQ. 8) THEN 

WRITE (11,11)  (TM(K,N) ,N=1. ORDER) 

WRITE (8, 83)  PS  (K) . (TM(K,n) ,N=1 , ORDER) 

FORMAT(I3, ' ] 1 ,8(13)) 

FORMAT  (812) 

ELSE 

WRITE (11 ,12)  (TM(K,N) ,N=1, ORDER) 

WRITE(8 , 93 )  PS  (K) , ( TM ( K , N ) ,N=1, ORDER) 

FORMAT ( 13 , 1 | 1 ,16(13)) 

FORMAT  (1613) 

ENDIF 

CONTINUE 

RETURN 

END 

SUBROUTINE  XMDAT ( TM , ORDER , THETA , NEXTST , PS , T , TD ) 

THIS  SUBROUTINE  FINDS  THE  NEXT  SYMBOL  STATE/S  BY  DETERMINING 
XMTR  MATRIX  CONSTRUCTED  BY  SUB.  XMSTA 

INTEGER  TM(64 , 64) ,PS(64) , TRDATA(256) , KK , K , THETA ( 64 ) , NN 
INTEGER  ORDER, T(256) ,TD(256) ,NEXTST(256) ,PRESTA(256) 

DO  5  1=1, ORDER 
T(PS (I ) )=I 
TD ( THETA ( I ) ) =1 
CONTINUE 
KK=0 

DO  7  K=0, ORDER- 1 
DO  8  M=0, ORDER- 1 
KK=KK+1 

NEXTST ( KK ) =TM ( T ( K ) ,TD(M) ) 

CONTINUE 

CONTINUE 

RETURN 

END 

SUBROUTINE  RCVR  (ORDER) 

THIS  SUBROUTINE  PREPARES  GENERAL  RECEIVER  MATRIX  AND 
SENDS  OUTPUTS  TO  THE  "FILE  RCDATA" . 

INTEGER  ORDER , MAT ( 64 ) , KL , RMAT (64,64) 

0PEN(UNIT=2,FILE=' RCDATA'  , STATUS= 1  OLD ' ) 

WRITE  (8,8l)0RDER 
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81  F0RMAT(///,T15, 'RECEIVER'S  DECODER  MATRIX  ANALYSIS 

*  FOR' ,13, ' -DPSK 1 , /, 

•'60( '-')  ,//,T5, 'GENERAL  RECEIVER  MATRIX1  ,/ ,T2,31(  1  -  1  ) ) 

DO  3  1=0, ORDER-1 
MAT ( 1+1 )=I 
3  CONTINUE 

DO  44  K=l, ORDER 
RMAT ( 1 , K) =MAT ( K ) 

44  CONTINUE 

DO  55  L=2, ORDER 
DO  66  M=1 ,L-1 

RMAT ( L , M ) =MAT ( ORDER-L+M+1 ) 

66  CONTINUE 

DO  77  N=L, ORDER 
RMAT ( L , N ) =MAT (N-L+l ) 

77  CONTINUE 
55  CONTINUE 

DO  88  NK=1, ORDER 
IF  ( ORDER. EQ. 3) THEN 

WRITE  ( 2 , 1 0 )  ( RMAT ( MK , KL ) , KL=1 , ORDER) 

WRITE  (8,15) (RMAT(NK,KL) ,KL=1, ORDER) 

10  FORMAT(8(I2) ) 

15  FORMAT ( 8X, 3 ( 12 ) ) 

ELSE 

WRITE  (2,11) (RMAT (NK,KL) , KL=1 , ORDER) 

WRITE  (8,16) (RMAT (NK , KL ) , KL=1 , ORDER) 

11  FORMAT (16 (12)) 

16  FORMAT (8X, 16(12)) 

END  IF 

88  CONTINUE 
CLOSE ( 2 ) 

RETURN 
END 

SUBROUTINE  RCVSTA  (PS , ORDER , THETA, RCVMAT) 

THIS  SUBROUTINE  PREPARES  SPECIFIC  RECEIVER  MATRIX  UNDER 
DEFINED  STATES. IT  READS  GENERAL  RECEIVER  MATRIX  FORM  FILE 
RCDATA  WHICH  HAD  BEEN  CREATED  BY  THE  SUBROUTINE  RCVDAT . 
OUTPUT  ARE  SENT  TO  THE  FILE  R. 


99 


199 


11 


12 


2 


13 
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INTEGER  TH(64,64) , ORDER, PS ( 64 ) , THETA (64) , RCVMAT (64, 64) 
OPEN (UNIT=2,FILE=' RCDATA' , STATUS3 ' OLD ' ) 

OPEN (UNIT=20 , FILE3 1 R '  , STATUS3 ' OLD 1 ) 

IF  ( ORDER. EQ- 8) THEN 

WRITE ( 8 , 99 ) (PS ( L ) , L=1 , ORDER ) , { THETA (I) ,1=1,  ORDER ) 
FORMAT(///,T5, 'RECEIVER  MATRIX  IS  1  , / , 31 ( 1  -  1 j , // . 

* 1  SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS ) 1 ,// ,4X, 8(1 
ELSE 

WRITE(8 . 199) (PS (L) ,L=1 , ORDER) 

FORMAT (///,T5, 'RECEIVER  MATRIX  IS ',/, 31 ('-'),// , 

*' SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS)' 

*,//,4X, 16(13)) 

ENDIF 

DO  1  1=1, ORDER 

IF  ( ORDER. EQ. 8) THEN 

READ (2,11)  (TH(I,N) ,N=1, ORDER) 

FORMAT (31 2) 

ELSE 

READ (2 ,12)  (TH(I,N) ,N=1, ORDER) 

FORMAT (1612) 

ENDIF 

DO  2  K=l, ORDER 

RCVMAT (I ,K)=THETA(TH(I,K)+1) 

CONTINUE 

IF  (ORDER. EQ.8)THEN 

WRITE (20, 13)  ( RCVMAT ( I , N ) ,N=1, ORDER) 

FORMAT( 812) 

WRITE (8, 33)  PS (I), (RCVMAT (I, N) ,N=1, ORDER) 

F0RMAT(I3, ' 1 ' ,8(13)) 
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ELSE 

WRITE(8 ,93)  PS(I) , (RCVMAT(I,N) ,N=1, ORDER) 

93  FORMAT  13, 1 1 ',16 (13)) 

WRITE (20, 14)  ( RCVMAT ( I ,  N ) ,N=1,0RDER) 

14  FORMAT ( 1612 ) 

END  IF 

1  CONTINUE 
CL0SE(20) 

RETURN 

END 

SUBROUTINE  RCVDAT (TR , PS , ORDER , NEXTST , RCVMAT ) 

THIS  SUBROUTINE  READS  SPECIFIC  RECEIVER  MATRIX  WHICH  HAD 
BEEN  CREATED  BY  SUBROUTINE  RCVSTA,THEN  SEARCHS  FOR  THE  RCVR 
OUTPUT . 


INTEGER  RCVMAT (64,64),PS(64),KK, ORDER 
INTEGER  TR ( 64 ) , NEXTST (256) 

0PEN(UNIT=20 , FILE* 1 R ' , STATUS* ' OLD ' ) 

DO  23  LL=1, ORDER 
IF  ( ORDER. EQ. 3) THEN 

READ  (20,24) (RCVMAT (LL, KM) ,KM=1, ORDER) 
24  FORMAT (8 (12) ) 

ELSE 

READ  ( 20 , 94 ) ( RCVMAT ( LL , KM ) , KM= 1 , ORDER ) 
94  FORMAT(16(l2)) 

END  IF 

23  CONTINUE 

DO  2  1=1, ORDER 
M=PS ( I ) 

TR(M)=I 

2  CONTINUE 

KK=0 

DO  3  M=0, ORDER- 1 
DO  4  K=0, ORDER -1 
KK=KK+1 

NEXTST (KK)=RCVMAT(TR(M) ,TR(K) ) 

4  CONTINUE 

3  CONTINUE 

RETURN 

END 


SUBROUTINE  CONVER (N , NEXTST , NN , ORDER , FL ) 


INTEGER  N,S,QU,SU,LLL,LLK(50) , MMM,NN , FL ,MM ,M , K(9 , 1024) 
INTEGER  MW,NT(256 , 4) , KZ , ORDER, NEXTST (2 56) ,ST,QQU,KD 


300  WRITE  (6 , *) 
WRITE  (6,*) 
WRITE  (6,*) 
WRITE  (6,*) 
WRITE  (6,*) 
WRITE  (6 .*) 
READ  (6,*)  S 
IF  (S.EQ.l) 
IF  (S.EQ.2) 


'PLEASE  ENTER  THE  SELECTION  CORRECTLY' 

'  ENTER  1  IF  YOU  NEED  BINARY  FORM  OF  SYMBOL  AND  ' 
'DATA  STATES  INCLUDING  SYMBOL  OUTPUTS' 

'  ENTER  2  IF  YOU  NEED  DECIMAL  FORM  OF  SYMBOL  AND 
'DATA  STATES  INCLUDING  SYMBOL  OUTPUTS' 

'  ENTER  3  IF  YOU  WANT  TO  TURN  BACK  MAIN  PROGRAM.' 


IF  (S.EQ.l)  GOTO  100 
IF  (S.EQ.2)  GOTO  200 
IF  (S.EQ.3)  GOTO  400 

IF  (S.NE.l .0R.S.NE.2.0R.S.NE.3)  GOTO  300 
IF  (FL.EQ.l)  THEN 
WRITE(8,6) 

FORMAT  (// , 8X, ' DATA  STATE ', 5X, 1  SYMBOL  STATE ', 5X, ' SYMBOL '  ) 
ELSE 

WRITE (8,21) 

FORMAT  (// , 8X, ' DATA  STATE ', 5X, ’ SYMBOL  STATE ', 5X, ' SYMBOL '  ) 
ENDIF 
KK=0 


nnnn 


DO  7  KZ=0, ORDER- 1 
DO  8  M=0, ORDER- 1 
KK=KK+1 

WRITE (8, 9)  KZ,M,  NEXTST (KK) 

9  FORMAT  (10X, 13 , 15X, 13 , 15X, 13) 

8  CONTINUE 
7  CONTINUE 
GOTO  300 

100  IF  (FL.EQ.l)  THEM 
WRITE (8, 47) 

47  FORMAT  (//.7X,'DATA  STATE ',10X, ' SYMBOL  STATE ',1 OX, 

* ' SYMBOL 1 , 72 ( 1  -  1 ) ) 

ELSE 

WRITE (3, 77) 

77  FORMAT  (//. 7X SYMBOL  STATE  1 , 10X, 1  DATA  STATE ',10X, 

*' SYMBOL' 72 ( ' - ‘ ) ) 

ENDIF 

I F ( ORDER . EQ . 8 ) THEN 

KD=3 

ELSE 

KD=4 

ENDIF 

DO  44  1=1 ,NN+1 
Q0U=i.Ea-'ST  (L ) 

DO  33  SR=KD ,1,-1 
NT ( L , SR ) =MOD ( QQU , 2 ) 

QQU=QQU/2 
33  CONTINUE 
44  CONTINUE 

DO  23  MW=0 ,NN 
QU=MW 

DO  22  LLL=N ,1,-1 
K(MW , LLL)=MOD(QU , 2 ) 

QU=QU/2 

22  CONTINUE 

IF  ( ORDER. EQ. 8) THEN 

WRITE (8,45)  MW, (K(MW,LK) ,LK=1 ,N) , (NT(MW+1 ,ST) ,ST=1 ,KD) 

45  FORMAT  (13 , 2X, 3(2X, 12 , 2X) , 3X,3 (2X, 12 , 2X) , 3X, 3 (2X, 12 , 2X) ) 

ELSE 

WRITE(8 , 55 )  MW, (K(MW,LK) ,LK=1,N) , (NT (MW+1 , ST) , ST=1 , KD) 

55  FORMAT  (13 , 2X , 4(2X, 12 , IX) , 3X, 4(2X, 12 , IX) , 3X, 4 (2X, 12 , IX) ) 

ENDIF 

WRITE  (8,46) 

46  FORMAT  ( 72 ( ' -  1 )  ) 

23  CONTINUE 

GOTO  300 
400  RETURN 
END 
C 

SUBROUTINE  ERROR ( FL , THETA , RCVMAT , PS , ORDER , TR , TM , T , NEXTST , TD ) 


THIS  SUBROUTINE  LOCATES  UNDETECTED  ERROR  COMBINATIONS. 


j", 

» 


INTEGER  QQU,QU . SOU,PS (64) , ORDER , RCVMAT (64,64) ,N,TR(64) ,TM(64,64) 
INTEGER  TD ( 2  5  6 ) ,  T  ( 2  5  6 ) , A , B , Y , KD , KN ( 1 ,8) , K J ( 1 , 8 ) , KB ( 1 , 8 ) , KY ( 1 , 8 ) 
INTEGER  FL, IS, NT (1,3) , KM ( 1 , 8 ) , NEXTST (256) , THETA (64) ,KH(1,8) 
OPEN(UNIT=l , FILE= 1  ERR '  , STATUS= 1  OLD  1 ) 

OPEN(UNIT=15 , FILE= 1 ERDATA ' , STATUS= 1  OLD 1 ) 

WRITE  (6  *)  1  IF  YOU  WANT  TO  HAVE  UNDETECTABLE  ERROR1 


WRITE (6,*)  'ANALYSES  ENTER 

READ  (6,*)  LD 

IF  ( ORDER. EQ. 8) THEN 

KD=3 

ELSE 

KD=4 

ENDIF 

IF (LD . EQ . 1 ) THEN 

WRITE  (6,*) 1  WHAT  KIND  OF 

WRITE  (6,*) 1  (1)  IN  BINARY 


, OTHERWISE  ENTER 


WHAT  KIND  OF  ANALYSES  OUTPUT  DO  YOU  WANT1 
(1)  IN  BINARY  FORM1 


*V  vj  i 


vv’VJ\°\  -  <■ 


»y»U  "H»  V 


13 

12 


333 

222 

444 

555 

45 

85 

55 

95 

10 


WRITE  (6,*)'  (2)  IN  DECIMAL  FORM1 

READ  (6,*)  IS 
WRITE (1,13) 

FORMAT/// //,72( '#’) ,// f7X, 1  CORRECT  RECEPTION  COMBINATIONS  IN  PRESE 
*NCE  OF  ERROR' ,//,72( '#' )) 

WRITE (1,12) 

FORMAT(// ,T5, 'INPUT' ,T17, 'PREVIOUS' ,T30, 'PRESENT' ,T43, 'PREVIOUS' , 
*T55 , 'PRESENT1 ,T63, 'OUTPUT' ,/,T15, 'TRANSMISSION' ,T23, 'TRANSMISSION' 
* , T42 , 'RECEPTION' ,T54, 'RECEPTION'  , / ,  79 ( ‘  - ' )) 

ELSE 

GOTO  111 
END  IF 
FL=2 

CALL  XMTR  (ORDER, FL) 

CALL  XMSTA  ( THETA, TM, PS , ORDER) 

CALL  XMDAT (TM , ORDER , THETA , NEXTST , PS , T , TD ) 

DO  1  1=0, ORDER-1 
DO  2  K=0, ORDER- 1 
A=TM(T(I ) ,TD(K) ) 

B=RCVMAT(TR(I ) ,TR(A) ) 

IF(IS.EQ.l) THEN 
QQU=A 

DO  333  JJ=KD ,1,-1 
NT ( 1 , JJ) =MOD (QQU , 2 ) 

QQU=QQU/2 

CONTINUE 

QU=K 

DO  222  LM=KD ,1,-1 
KM( 1 , LM)=MOD (QU , 2 ) 

QU=Q U/2 
CONTINUE 
SQU=I 

DO  444  LN=KD ,1,-1 
KN ( 1 , LN ) =MOD ( SQU , 2 ) 

SQU=SQU/2 

CONTINUE 

SQU=B 

DO  555  LN=KD ,1,-1 
KB ( 1 , LN ) =MOD ( SQU , 2 ) 

SQU=SQU/2 

CONTINUE 

IF  (ORDEP  SO. 8) THEN 

WRITE (1 ,-.j) (KM(1 ,  L) ,L=1.KD), (KN(1,KK) ,KK=1,KD) , (NT ( 1 , IT) , IT=1 , KD ) 
*, (KN(l,KK),KK=l,KD), (NT ( 1 , IT) , IT=1 , KD) , (KB(1,IJ) ,IJ=1.KD) 

FORMAT (T3, 3 (13) ,T15 ,3(13) ,T28,3(I3) ,T42,3(I3) ,T54,3(I3) ,T65,3(I3) ) 
WRITE (15, 85) (KN(1 ,KK) ,KK=1 , KD) , (NT(1 , IT) ,IT=1 ,KD^ 

*,(KB(1IJ)  IJ=1,KD) 

FORMAT(9(I2) ) 

ELSE 

WRITE (1, 55 ) (KM(1,L) ,L=1 .KD) , (KN( 1 , KK) , KK=1 , KD)  (NT(1 , IT) ,IT=1 . KD 
*, (KN( 1 ,KK) , KK=1 , KD)  , (NT ( 1 , IT) ,IT=1,KD) , (KB (l , I J ) , I J=1 . KD ) 

FORMAT (T1 ,4(12) , T15 ,4(12) ,T28 ,4(12) ,T42 ,4(12) ,T54 , 4(12 ) ,T65 .4  If 
WRITE (15, 95) (KN(1 ,KK) ,KK=1 ,KD) , (NT( 1 , IT) , IT=1 , KD) 

*,(KB(1  IJ),IJ=1,KD) 

FORMAT/ 12 (12)) 

END  IF 
ELSE 

WRITE (1,10)  K , I , A , I , A , B 

FORMAT (T 10 ,I3,T25,I3,T35,I3,T48,I3,T58,I3,T71,I3) 

END  IF 

DO  3  L=0, ORDER- 1 

DO  4  J=0, ORDER-1 

I F ( L . EQ . I . AND . J . EQ . A ) THEN 

GOTO  100 

ELSE 

Y=RCVMAT(TR(L ) ,TR(J) ) 

IF(Y . EQ . B)THEN 
IF( IS .EQ . 1 )THEN 
SQU=L 

DO  666  LN=KD ,1,-1 

95 
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KH(1 ,LN)=MOD(SQU, 2) 

SQU=SQU/2 

CONTINUE 

SQU=J 

DO  777  LN=KD ,1,-1 
K J ( 1 , LN ) =MOD ( SQU , 2 ) 

SQU=SQU/2 

CONTINUE 

SQU=Y 

DO  888  LN=KD ,1,-1 
KY ( 1 , LN ) =MOD ( SQU , 2 ) 

SQU-SQU/ 2 
CONTINUE 

IF(ORDER.EQ.8)THEN 

WRITE ( 1 ,48) (KH( 1 ,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=1,KD) 

* , (KY(1 , IT) ,IT=1 ,KD) 

FORMAT  (T42 ,3(13) ,T54 , 3 ( 13) ,T65 , 3 (13 ) ) 
WRITE(15,999)(KH(1,KK) ,KK=1,KD) , (KJ(1,KK) ,KK=I,KD) 
* » (KY ( 1 , IT) , IT=1 , KD ) 

FORMAT  (9(12)) 

ELSE 

WRITE (1,49)(KH(1,KK) ,  KK=1 , KD) , (KJ(1 ,KK) , KK=1 ,  KD) 

*, (KY(1,IT) ,IT=1.KD) 

FORMAT  (T42 ,4(12) ,T54,4(I2) ,T65,4(I2)) 

WRITE ( 15 , 949) (KH(1 , KK) , KK=1 , KD) ,  (KJ(1 , KK) , KK=1  ,KD) 
* , ( KY ( 1 , IT ) , IT=1 , KD ) 

FORMAT  (12(12)) 

END  IF 
ELSE 

WRITE (1 ,20)  L,J,Y 
F0RMAT(T48, 13 ,T58, 13 ,T71 , 13) 

END  IF 
END  IF 
END  IF 
CONTINUE 
CONTINUE 
WRITE  (1,46) 

FORMAT  (79( ' -  1 ) ) 

CONTINUE 
CONTINUE 
CLOSE ( 1 ) 

CLOSE (15) 

RETURN 

END 


v  v  --  .*  •->  ‘j 


APPENDIX  B 

COMPUTER  GENERATED  ENCODER  AND  DECODER  OUTPUTS 


1.  8-DPSK  TRANSMITTER 


TRANSMITTER'S  ENCODER  MATRIX  TABLE  FOR  8-DPSK 
GENERAL  TRANSMITTER  MATRIX 


ABCDEFGH 

BCDEFGHA 

CDEFGHAB 

DEFGHABC 

EFGHABCD 

FGHABCDE 

GHABCDEF 

HABCDEFG 


TRANSMITTER  MATRIX 


SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS) 
/  01326754 


75132046 

51320467 


13204675 

32046751 


20467513 

04675132 


46751320 

67513204 


SYMBOL  STATE 
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SYMBOL 

0 
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0 
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1 
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5 

0 

0 

0 

1 

0 

1 

0 

1 

1 

6 

0 

0 

0 

1 

1 

0 

1 

0 

1 

7 

0 

0 

0 

1 

1 

1 

0 

0 

1 

8 

0 

0 

1 

0 

0 

0 

0 

0 

1 

9 

0 

0 

1 

0 

0 

1 

0 

1 

1 

10 

0 

0 

1 

0 

1 

0 

0 

0 

0 
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0 

0 
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0 

1 

1 

0 

1 

0 

12 
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0 

1 

1 

0 

0 

1 

0 

1 

2.  8-DPSK  RECEIVER 


RECEIVER'S  DECODER  MATRIX  TABLE  FOR  8-DPSK 


GENERAL  RECEIVER  MATRIX 


01234567 

70123456 

67012345 

56701234 

45670123 

34567012 

23456701 

12345670 


RECEIVER  MATRIX  IS 


SYMBOL  STATE  (COLUMNS)  /  SYMBOL  STATE  (ROWS) 
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5 
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2 

0 
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6 
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1 

3 
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7 

5 

4 

0 

7 
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1 
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6 

7 

5 

4 

5 

4 
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3 

2 

6 

7 

5 

1 

5 

4 

0 

1 

3 

2 

6 

7 

3 

7 

5 

4 

0 

1 

3 

2 

6 

2 

6 

7 

5 

4 

0 

1 

3 

2 

0 

2 

6 

7 

5 

4 

0 

1 

3 

4 

3 

2 

6 

7 

5 

4 

0 

1 

6 

1 

3 

2 

6 

7 

5 

4 

0 

SYMBOL 

STATE 

SYMBOL  STATE 

DATA  STATE 

51 

1 

1 

0 

0 

1 

1 

1 

1 

0 

52 

1 

1 

0 

1 

0 

0 

1 

0 

0 

53 

1 

1 

0 

1 

0 

1 

0 

1 

1 

54 

1 

1 

0 

1 

1 

0 

0 

0 

0 

55 

1 

1 

0 

1 

1 

1 

0 

0 

1 

56 

1 

1 

1 

0 

0 

0 

1 

1 
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1 

1 

0 

0 

1 

0 

1 

1 

53 

1 

1 
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1 

1 

1 

0 

1 

1 

0 

1 

0 

60 

1 
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0 
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0 

1 

62 

1 

1 

1 

1 

1 

0 

1 

0 

0 

63 

1 

1 

1 

1 

1 

1 

0 

0 

0 

TRANSMITTER 

INFORMATION 

FOR  THE 

GIVEN 

RECEIVER 

GENERAL  TRANSMITTER  MATRIX 


ABCDEFGH 

BCDEFGHA 

CDEFGHAB 

DEFGHABC 

EFGHABCD 

FGHABCDE 

GHABCDEF 

HABCDEFG 


TRANSMITTER  MATRIX 


SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS) 
01326754 

7  75132046 

5  51320467 

1  13204675 

3  32046751 

2  20467513 

0  04675132 

4  46751320 

6  67513204 


3.  16-DPSK  TRANSMITTER 


TRANSMITTER'S  ENCODER  MATRIX  TABLE  FOR  8-DPSK 


GENERAL  TRANSMITTER  MATRIX 


ABCDEFGHI JKLMNOP 
BCDEFGHI JKLMNOPA 
CDEFGHI JKLMNOPAB 
DEFGHI JKLMNOPABC 
EFGHIJKLMNOPABCD 
FGHI JKLMNOPABCDE 
GHI JKLMNOPABCDEF 
HI JKLMNOPAB CDEFG 
I JKLMNOPABCDEFGH 
JKLMNOPABCDE  FGHI 
KLMNOPAB  CDE  FGHI J 
LMNOPABCDEFGHIJK 
MNOPABCDEFGHI JKL 
NOPABCDEFGHI JKLM 
OPABCDE  F  GHI JKLMN 
PABCDEFGHI JKLMNO 


TRANSMITTER  MATRIX 


SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS) 
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13 
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0 

2 
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8 
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14 
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13 

9 

11 

3 

1 

5 

7 

6 

4 

2 

2 

10 

8 

12 
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15 
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5 
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6 

4 

0 
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10 
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14 

15 

13 

9 

11 

3 

1 

5 

7 

6 

4 

0 

2 

8 

8 

12 

14 

15 

13 

9 

11 

3 

1 

5 
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6 

4 

0 

2 

10 

12 

12 

14 

15 

13 

9 

11 

3 

1 

5 
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4 

0 

2 

10 

8 

14 

14 

15 

13 

9 

11 

3 

1 

5 

7 

6 

4 
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SYMBOL  STATE 


DATA  STATE 


SYMBOL 
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4.  16-DPSK  RECEIVER 


RECEIVER'S  DECODER  MATRIX  TABLE  FOR  16-DPSK 


GENERAL  RECEIVER  MATRIX 


012345678  9101112131415 
15  012345678  91011121314 
1415  012345678  910111213 
131415  012345678  9101112 
12131415  012345678  91011 
1112131415  012345678  910 
101112131415  0123456789 
9101112131415  012345678 
8  9101112131415  01234567 
7  8  9101112131415  0123456 
678  9101112131415  012345 
5678  9101112131415  01234 
45678  9101112131415  0123 
345678  9101112131415  012 
2345678  9101112131415  0  1 
12345678  9101112131415  0 


RECEIVER  MATRIX  IS 


SYMBOL  STATE  (COLUMNS)  /  SYMBOL  STATE  (ROWS) 

15  13  9  11  3  1  5  7  6  4  0  2  10  8  12  14 


15  0  1  3  2  6  7  5  4  12  13  15  14  10  11  9  8 

13  8  0  1  3  2  6  7  5  4  12  13  15  14  10  11  9 

9  9  8  0  1  3  2  6  7  5  4  12  13  15  14  10  11 

11  11  9  8  0  1  3  2  6  7  5  4  12  13  15  14  10 

3  10  11  9  8  0  1  3  2  6  7  5  4  12  13  15  14 

1  14  10  11  9  8  0  1  3  2  6  7  5  4  12  13  15 

5  15  14  10  11  9  8  0  1  3  2  6  7  5  4  12  13 

7  13  15  14  10  11  9  8  0  1  3  2  6  7  5  4  12 

6  12  13  15  14  10  11  9  8  0  1  3  2  6  7  5  4 

4  4  12  13  15  14  10  11  9  8  0  1  3  2  6  7  5 

0  5  4  12  13  15  14  10  11  9  8  0  1  3  2  6  7 

2  7  5  4  12  13  15  14  10  11  9  8  0  1  3  2  6 

10  6  7  5  4  12  13  15  14  10  11  9  8  0  1  3  2 

8  2  6  7  5  4  12  13  15  14  10  11  9  8  0  1  3 

12  3  2  6  7  5  4  12  13  15  14  10  11  9  8  0  1 

14  1  3  2  6  7  5  4  12  13  15  14  10  11  9  8  0 


SYMBOL  STATE 


SYMBOL  STATE 


DATA  STATE 


TRANSMITTER  INFORMATION  FOR  THE  GIVEN  RECEIVER 


GENERAL  TRANSMITTER  MATRIX 


ABCDEFGHI JKLMNOP 
BCDEFGHI JKLMNOPA 
CDEFC-H I  JKLMNOP  AB 
DEFGHI JKLMNOPABC 
EFGHI JKLMNOPABCD 
FGHI JKLMNOPABCDE 
GHI JKLMNOPABCDEF 
HI JKLMNOPABCDE FG 
I JKLMNOPABCDEFGH 
JKLMMOPABCDEFGHI 
KLMNOPABCDEFGHI J 
LMNCPABCDEFGHI JK 
MNOPABCDEFGHI JKL 
NOPABCDEFGHI JKLM 
OPABCDEFGHI JKLMN 
PABCDEFGHI JKLMNO 


TRANSMITTER  MATRIX 


SYMBOL  STATE  (COLUMNS)  /  DATA  STATE  (ROWS) 
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APPENDIX  C 

COMPUTER  GENERATED  OUTPUTS  FOR  DBDR'S  ERROR  PATTERN 

ANALYSES 

The  program  outputs  given  below  are  a  partial  sample  of  the  main  outputs. 
Since  undetectable  error  patterns  are  similar,  in  the  next  sections,  1  8  period  of  the 
pattern  have  been  included.  Actually  526  combinations  of  error  patterns  are  detected 
by  the  program  for  the  S-DPSK  case.  For  the  16-DPSK  case  there  are  4096 
undetectable  error  pattern  combinations. 

1.  8-DPSK  RECEIVER 

CORRECT  RECEPTION  COMBINATIONS  IN  PRESENCE  OF  ERROR 


INPUT  PREVIOUS  PRESENT  PREVIOUS  PRESENT  OUTPUT 

TRANSMISSION  TRANSMISSION  RECEPTION  RECEPTION 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

1 

1 

0 

1 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

1 

1 

0 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

I 

1 

0 

1 

0 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

0 

1 

1 

0 

1 

0 

0 

1 

0 

0 

1 

1 

I 

0 

1 

1 

1 

0 

0 

I 

1 

1 

1 

1 

0 

1 

0 

0 

1 

0 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

I 

0 

0 

1 

0 

1 

1 

0 

0 

i 

0 

0 

1 

1 

1 

0 

0 

0 

1 

0 

1 

0 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

0 

1 

0 

1 

1 

0 

0 

0 

1 

0 

I 

0 

1 

1 

1 

0 

1 

1 

0 

1 

0 

0 

I 

1 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

0 

1 

0 

0 

1 

I 

0 

1 

0 

1 

0 

0 

0 

1 

1 

0 

1 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

1 

1 

1 

0 

1 

I 

1 

0 

1 

0 

1 

1 

0 

1 

1 

1 

1 

0 

1 

0 

1 

0 

1 

1 

1 

1 

1 

0 

0 

1 

0 

1 

1 

1  0  0  0  0  0 


0  10 


0  0  0 
0  0  1 
0  10 
0  1  1 
10  0 
1  0  1 
110 
111 


0  10 
1  0  1 
0  1  1 
0  0  1 
0  0  0 
1  1  1 
10  0 
1  1  0 


10  0 
1  0  0 
1  0  0 
1  0  0 
10  0 
10  0 
10  0 
10  0 


■■■ 


2.  16-DPSK  RECEIVER 


######################################################################## 
CORRECT  RECEPTION  COMBINATIONS  IN  PRECENSE  OF  ERROR 
St###################################################################### 


INPUT 


PREVIOUS  PRESENT 
TRANSMISSION  TRANSMISSION 


0  0  0  0 


0  0  0  0 


0  0  0  0 


PREVIOUS 

RECEPTION 


0  0  0  0 
0  0  0  1 
0  0  10 
0  0  11 
0  10  0 
0  10  1 
0  110 
0  111 
10  0  0 
10  0  1 
10  10 
10  11 
110  0 
110  1 
1110 
1111 


PRESENT 

RECEPTION 


0  0  0  0 
0  0  0  1 
0  0  10 
0  0  11 
0  10  0 
0  10  1 
0  110 
0  111 
10  0  0 
10  0  1 
10  10 
10  11 
110  0 
110  1 
1110 
1111 


OUTPUT 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


0  0  0  1 


0  0  0  0 


0  0  10 


0  0  0  0 
0  0  0  1 
0  0  10 
0  0  11 
0  10  0 
0  10  1 
0  110 
0  111 
10  0  0 
10  0  1 
1  0  1  ‘ 
10  1 
1  1  0 
110 
1  1  1 
111 


0  0  10 

10  10 
0  0  0  1 
0  0  0  0 
0  111 
0  10  0 
0  110 
110  0 
10  11 


0  0  0  1 


0  0  0  1 
0  0  0  1 
0  0  0  1 
0  0  0  1 
0  0  0  1 
0  0  0  1 
0  0  0  1 
0  0  0  1 


§ 


£ 


0  1 
0  110 
0  111 
10  0  0 
10  0  1 
1  0 


110  0 
10  0  0 
10  11 
0  110 


0  10  1 
0  10  1 
0  10  1 
0  10  1 


10  11 
110  0 


10  0  1 
10  0  0 


1110 


10  0  0 
10  0  0 


10  0  0 


0  0  1 
0  0  1 
0  10 
0  10 
Oil 
Oil 
10  0 
10  0  1 
10  10 
10  11 
110  0 
110 
111 
111 


1110 


0  0  0  0 


0  0  0  1 


0  0  0 
0  0  0  1 
0  0  10 
0  0  11 
0  10  0 
0  10  1 
0  110 
0  111 
10  0  0 
10  0  1 
10  10 
10  11 
110  0 
110  1 
1110 
1111 


0  0  10 
0  0  0  1 
10  10 
0  111 


1111 
0  10  1 
1110 
0  0  11 
110  1 
10  11 
10  0  1 
0  110 
10  0  0 
0  111 
110  0 
0  10  0 
10  10 
0  0  0  0 
0  0  10 


110  1 
110  1 
110  1 
110  1 


1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 

1110 


1111 


0  0  0  0 


0  0  11 


0  0  0  0 
0  0  0  1 
0  0  10 
0  0  11 
0  10  0 
0  10  1 
0  110 
0  111 
10  0  0 
10  0  1 


0  0  11 
1110 
0  0  0  1 
110  0 
10  11 
1111 
10  0  1 
110  1 
0  111 
10  10 


1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 
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